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ABSTRACT 


One  of  the  main  problems  connected  with  the  analytical 
evaluation  of  primary  loads  and  transient  control  effects  on 
helicopter  rotors  consists  in  the  assumption  of  the  local 
induced  velocities  over  the  rotor  disc. 

This  problem  has  attracted  the  attention  of  many  engineers 
all  over  the  world  but,  up  to  now,  no  satisfactory  engineering 
solution  w as  found  and  though  the  mean  value  of  the  induced 
velocity,  in  hoverin,j  and  forward  flight  conditions,  is  well  # 
known,  thus  permitting  a  fairly  good  evaluation  of  the  helicojli  ' 
ter  general  performance  characteristics,  lack  of  knowledge  of  th« 
actual  vortex  structure,  which  gives  rise  to  the  local  blade 
induced  velocities,  explains  the  poor  agreement  between  analyti¬ 
cal  and  experimental  results,  mainly  in  blade  stress  problems 
and  those  connected  with  the  transient  effects  in  rotor  control 
and  stability  studies. 

The  solution  of  this  problem  by  analytical  methods  alone 
leads  to  practically  insoluble  difficulties  and  for  this  reason 
an  experimental  method  consisting  in  visualizing  the  vortices 
by  smoke  emission  on  model  rotors , operated  in  a  wind  tunnel, was 
used  to  help  to  build-up  a  reasonable  and  intentionally  simpli¬ 
fied  (for  mathematical  accessibility)  canevas  for  the  rotor 
vortex  system. 

As  a  result  of  this  study  engineering  methods  are  proposed 
for  the  analytical  evaluation,  with  a  satisfactory  degree  of 
approximation,  of  the  local  instantaneous  induced  velocities  over 
a  rotor  disc,  in  various  helicopter  flight  conditions,  as  a  func¬ 
tion  of  the  rotor  characteristics. 

Some  of  these  methods  are  immediately  applicable,  others 
need  numerical  function  tables  which  remain  to  be  established. 
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LIST  OP  SYMBOLS 


R 


c 

a 


Rotor  radius,  measured  from  center  of 
rotation 

distance  of  a  blade  section  to  center  of 
rotation 


(m) 


R 


r  - 


R 


•V ,  jk. 

*  •  • 

distance  of  inner  end  of  blade  to  center  of  '*■'  w  *• 
rotation  '  #  *(»?'•* 

*  ,  • 


*• 

•t 

i. 


dummy  variable  of  blade  irlemenfc 

•  ,  < 

blade  section  chord 

% 

distance  of  flapping  hinge  to.  center  of 
rotation 


• 


b 

k 

6c 

0o6a  • 

i 

V  ..  *' 


constant  coefficient 
number  of  blade* 
rotor  disc  aern  *. 


•  •  r 


.  »  . 


&  a* 


V*  * 


*  .  •(») 

.  *.  ■* 

•  • 

**  -.<m) 

•  *•  » 

'  * 

•»  ..«* 

*  * 

•*!»  (sq.m.) 


blade  section  pitch  angle 


"M 


cyclic  pitch  angles  ..  ..»< 

•  „  *  c.-  .  :sv-  •  .**  4  •, 

bladq  tv.’ist  ahf*l 

0  -A*  * 


©c*» 

•  A  ^ 


blade  section  angle  of  attack'o-c>*.-' 

.•••*«  f  «*  .  V. : 

blade  .azimuth  angle  ’ 


rotor-shaft  tilting  angle  (positive  nose  do;#n 
unless  otherwise  indicated  in  text. 

See  paragr.3 .2M 

angle  of  rotation  of  airfoil  section 
blade  flapping  angle 
flapping  Fourier  series  coefficients 
rotor  solidity 
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period  of  rotation 


blade  circulation 


blade  Lock's  constant 
blade  aerodynamic  lift 
*  rotor  thrust 


,  ^  Section  lift  ca«f f ioiertt  '*• 

•*  '  *  !  '  >  F  \ 

•  #  rotor  thrust  coef£ icleitt'-  (  =  ^Cr) 

**  Tm  rotcr  drag  *.  ...  *  4j’  ‘  •  r~- 
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V.  induce^  vnlocity  on  trie  olade  element  .  .  . 

1  ,  -  •-•  t4  ^  •  *  •  •  «  A.  .> 

Vf.Vi  ind«c#d^flo»  '(bVneath  rotor  discj  .  *  *  • 

^  •  •  "J4-*/  "  mf  JM  <>*  , 
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£  I  blade  flexural  ri<Ji<i!*y  c  *  *  '*  .  * 

kl  •  .*  ..•  •  -  -  .  ’  ■  ''',  ^  \, 

y  blade  natural  t  rcduen^jr  '•  •  **•  '  •..  '  *' 

'  -  *  ...  *  *■•  ..*•<*  ’■  . 


i  *•  3«  blaie  shade  parcmetvr 

*  • 


'sn*  fu&r'dr. 

*  •*,  #  *  1  z  dT  <; 

•  J  a  ,  -  .  . 


Moment  af  inertia  of  blade  rbOTH; 
fleppin^  hinge  **;V^  •  *  .- 


All  other  particular  symbols  used  are  explained 
in  the  text. 
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EXPERIMENTAL  AND  THEORETICAL  STUDY 
OF  LOCAL  INDUCED  VELOCITIES  OVER  A  ROTOR 
DISC  FOR  ANALYTICAL  EVALUATION  OF  PRIMARY 
LOADS  ON  HELICOPTER  ROTOR  BLADES 


i 
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The  purpose  of  this  stud/  is  to  make  available  to  engineers 
a  method  of  evaluating,  with  a  satlsfaotory  degree  of ,approxima- . 
tlon ,  the  aerodynamic  loads  distributed  over  the  blades  of  a  heli¬ 
copter,  rotor.  ...  l-1  "« 


t*  K->*  *  **  ‘  *  *  *  ‘  ,*'•  •  ?  ” 

1  .-cs  Two  aaln  off  sets  were  introduced  c 

Y  ’  ,  •  ••  :  a  '  „  ^.0c;5 


'  5  V 

„  .■  *-0  c'  r  j  t'/  v-.  J  v"',» 


a  m  ,‘v .  $y 

4 *'• 


-•"velocity  induced  by  the  free  vortices.' 


-  non  stationary  regime. 


°  *-v'  ‘-4  ^ 


This  doounent  Is  divided  Into  two  parts  t 
-  an  experimental  part. 


-  an  analytical  part. 


OO  ^  Cr  w  v-  -• 

c>-  •  <rJc ^ 


In  the  experimental'  part,  an  attempt*  war;hede .  to-  situate 
In  space  the  free  vortices  issuing  from  ihma  roto^f or  different 
helicopter  flight  configurations.  Indeed,  , in jtbe^dase/of  fixed 
wings,  the  location  of  the  free  vort lees  Is  known,  whereas  It  is 
difficult  to  define  theoretically  the  location' of  suohvortices 
beneath  the  rotor  disc  in  the  case  of  heilooptor.rof orar  it 
was  therefore  decided  to  resort  to  experiment,- which; consisted 
In  Visualizing  the  vortices  by  means  of  ancle.  On'  model,  rotors 
operated  in  a  wind  tunnel. 


In  the  analytical  part,  the  proposed  method  "Is “set  out  and 
illustrated  with  hUBerioal  examples.  . 


■  v  -v  ,, 


.The  method  is  an  approximate  one  CftlyV  'Wherever  possible 
conventional  calculation  metheds  have  been  use  3 


WIND-TUNNEL  VISUALIZATION  TEST  PROCEDURE 


-  Description  of  method 


The  visual izat Ions  were  made  in  order  to  locate  the  free 
vortices  streaming  from  the  blades. 
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A  double  visualization  was  performed  r 

-  First,  the  smoke  Jet  was  emitted  upstream  and  above  the 
swept  disc.  These  locations  were  chosen  so  that  the  Jet 
touched  the  tip  of  the  blade  for  a  given  asimuthal  posltlo 
(see  photographs  Fig.  3  through  10).  .  This  provided  sec¬ 
tion  views  of  the  vortex. 

-  Second,  the  smoke  was  emitted  either  at  the  two  blade 
tips  or  at  the  root,  thereby  showing  up  the  vortices  in 
continuous  fashion  (see  photographs  Fig. 11  through  16). 

The  lift,  the  flapping  and  the  drag  were  recorded  in  each 
case. 

For  the  dynamic  regimes,  the  input  force  was  also  recorded, 
but  on  the  other  hand  the  drag  measurement  was  suppressed. 

The  tests  were  performed  at  the  wind-tunnel  of  the  INSTITUT 
DE  MECANIQUE  DES  FLUIDES  at  MARSEILLES,  under  the  direction 
of  Professor  VALENSI. 

2  -  Description  of  Rotor 

1.21  -  Definition  of  Blades 


1.211  -  Basic  Characteristics  (See  tables  1.1,  1.2  and  1.3) 

In  order  to  facilitate  data  processing  and  interpre¬ 
tation  the  blades  were  of  rectangular  plan  form  (cons¬ 
tant  chord),  without  twist.  The  airfoil  adopted  was 
the  NACA  0018.  Rotor  radius  was  R  «  0.75m(2.46  ft). 

From  the  structural  standpoint,  the  blades  comprised 
a  massive  leading  edge  spar  (see  Figure  47)  and  two 
shells  in  thick  plate. 

It  was  sought  to  make  the  blades  rigid  in  bending  and 
torsion  (see  table  1.3)  to  avoid  introducing  the  bon¬ 
ding  and  pitch  distorsions  in  the  calculations. 

To  avoid  loads  or.  the  controls  and  a  parasite  cyclic 
pitch  tfe  ,  the  C.G.  line  was  brought  as  close  as  1 

possible  to  the  feathering  axis.  Blade  No.  5  deviates 
most  from  this  condition  (see  Table  1.1). 

1.212  -  Specific  Test  Requirements 

To  permit  photographing  to  take  place,  it  was  necessarj 
for  the  rotor  to  be  able  to  rotate  in  both  directions. 
This  was  possible  by  reversing  the  blades,  due  to 
the  fact  that  they  were  symmetrical  (NACA  0018  airfoil, 
with'  no  twist) . 

/. 
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P°r  the  smoke  emissions  from  the  blade  tips,  a  stain¬ 
less  steel  tube  was  incorporated  in  the  blade  (Pig. 47) 


The  blades  were  balanced  statically  and  dynamically 


_Ratl»_and  of  Look's  Constant 


Tha^dif farantyaeta  ^  cfnditlon#  achieve*  with  the 
ef'fijides  t  .  %  ,  ’  ,  e*  ■  %  * 

Oi£Up  1  Xb=id*|  1,  V:,3*jVioIsinyir*  -  0*0*7*11 

■  ■  Vs  ^  ■ .« 

i,X  iJ  -I  >  ./,f  v, 

, ^  -  ■  i-v  r^i-  T  r’4'  ^  .  'f  1** 

6  c£?  .  life  */&s**’  fei,,  j  *  "A,  &  «Ctt  ,V  #  «  **  1  ‘ 


■§lad  e  p^^daftfga'ed 


IB- subh/a  wajKfhat  thalr  natural* 
torsi onar.pnd^, bend! ag  freQuehoies/were  widely  d’istinet. 
.  i^^er  t^dhviaete  b e ndlng/tors  ion vecupl  mgs  .0,. tv  ^ 

*■  v? 

^•r®  alio  oonp&red  id  tho  fretjuenoie 
aerodynamic,  forces  (up  to  lOcu)  so  that  there 
fehouia  toe  ho  reaonanee  at  the  chosen  rotational  anted  ; 

.s’  ?a>l?in8  i/lnto oA^e'ouftt  thefeot  thatthe,  blades  were  ,i‘f 

v-.(  beawy  4ndr  that  it  ,was  also  necessary  to  insure  e  •  k  . 

•  ,  ,re_ftixire»entni  good  teyaienal  r  if  Wit#  in*  eonitanoy^ 
;o-  -0.»r-the  problen  was  obviously  a  tricky  one.  :<&'•:  . 

fi -r^Sfikpe  Parameters  „cm\v  £ l- JVr <’£  ♦  •’■ D " X?. *&  •  • 


Jjaechlade  was  defined'  bySls.lujf^pn^iSBet'e^'s^f  the 
■oiauert  form  (table  *:£?.*%  .^v*.  ‘r.v  ,.* 
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Por  a  rectangular  blade  t 
.  Go  ^  1 

G1  5s?  1/2 

•  »i 

tfp  1/5"* 


if  B 


•  ’*  ,o  * 

1*22  -  ftotor*  •"Hub*  -  Control^  ^  ^  f4  ‘ . 


'  4* 


* •  «p*t? ?* 


•  .  ,  i 

1.221  -  The  rotor  mat  flOuntad  OD  the  ecd  Of  tho  balance  arm 
(photographs  Fig.  1  and  2)  provided  *tth  oyolic  pit  oh 
(  and  t oolltotiv*  pitch  (©«.)  oontrols  oapable 

*.  of  reaott  control .  it  was  also  posslbl#  to  tilt  the 

rotor  in  pitoh-  (<*  )  by’rsaote  .control.  * 


■The  rotor  was  drift©'  ‘by  an  asynchronous  motor. 


The  shaft  passed  through  the  balance  arm'* 


’  ».  • 


•  V?  *  *■'  ’«»  •  *  <  A  4i^  «>»  .  , 

Snoice  was  de  live  rod,  t'o  the  blados.  through  ’the 
,6OOn.tor  of  the/*hubi  7* •  „•  •  .. „ 

*  **>  '1?X.  !  <•  J‘  .  ■*“ v  •  c:'  .  tf*'  * 


■’  l«222  » •  Balancing 


«*  ■  *  •  .  t  •  _  «-  ^  rv  V 

All" rotating'  control  'm^aseo -  (levers.,  and  swaahplate) 

.  had  to  bo •  Ve cy "’Oa fjfuJ Ilk. ha i ane «d  statically  and  dyna- 

■  mloall j. . % .'W;,-, i" r.  ‘  • 

•  **  f  ;  •*  '•  ^  c  ‘J  .<»  t.. .>*  a.  *  ^  ^  *5>»  *  •  • 

,  ...  ^  .•  n  <*>  ;/  ,  «-*>,. 

Tho  natural  frYquenei’o'i  rp,r  the  d rive  shaft  and  of  the 
tunhfl  bala'aca jraTe  compared  to  thoao  of.  the  blades 
(Fig. 55).  *  "  ^  v  •* 


f?»  wVi^y; 


1.5  last  .Program 
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•'  i* 
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_  AotOltfapijO  a/i<i  rtoordlr^ii-stro^iade  for  groups  A  and  B, 


Vilms  And  feoordthgs* sore  Vide  for  group  c. 
"  firoupa  a  and,  B . 

1.511*  Flight  Parameters 


It  wa^  desirable  to  ’cover  different  flight  conditions  in 
order  to  be  able  to  locate  the  vortices  for  each  different 
configuration.  To  this  end,  the  following  were  made  to 
vary  i  / 

•  /  • 


Table  1. 
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Torsional  Rigidity 
Flexural  Rigidity 


Blade  n° 


Torsion  0  0 


for  Fw  =  20  Kg. 


Bending 


(Kg. m2) 


20/100  20/100  15/100  20/100 


87,2 


13/100 


•  m  <m  f  *w *  «  ,  * 

.*  • 

:  -  ••  4' 

«?V  •>  **  " 

*'•’  f;  .<tj, 

.  V  -‘i  *  if  .  -  ' . 

4i  wA  •«  * >v-  •'  ' 

•  W  , 

■>  V 
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-  the  collective  pitch 


-  the  cyclic  pitch  Q.  Q. 

•  •  V  •  • 

0  4  6C  ♦  ®  ^  +  0a  -S«n  yf 

•i  .  • 

*•  V  the  inclination  of  the  j'gtor  shafg 

*  .*•  \  V.  •  • 

•  %  «  '**• 

-  the,  sol  idi4f.  .  .  .  .  »• 

■:»  »«  .  *  . 
•-‘LOdk's  C6fcstjifl%  0  * 

»  *  •  t  • 

Xhe  v*iue»  of  parameters  and  their  variations 

.  veri  determined  £©  as  ift  fulfil  both  mechanical  end  pho¬ 
tographic  Requirements . 


This  vAe  important  because/  from  the  standpoints  of  pro* 
cessing  and, visualisation, . any  tiro  values  of  a  given  para 
.■step  had  to  differ  sufficiently  from  each  other-eo  that 
‘  ”&  priori"  the  effects  should  also  be  different  and  the 
vortices  distinct  from  each  other* 

v  •  ,  •  ‘  i  .  •. 


Blade  No,  4  (o*  0.118  tn  or  i0*?9  ft)  was  chosen  as  the 
basic  blade  for  the  teats  as  it  best  satisfied  the  rigi¬ 
dity*  C*0.  and  weight  requirements* 


*.*  »  .  .  •  •• 

•Test  Phases 

•  ...  . 

Tables  1*4  and  1.5  (pages  A.  10 through  A. 36)  indicate  the 
test  progress  for  these  two  groups. 

s  *  * 

•  "  h  ,  « 

Each  individual  case  was  itself  divided  into  several 
groups  of  photographs.  It  is  proposed  to  examine  one 
epeolfie  esse  in  detail. 


External  Smoke  Emission  (Photographs  3  through  10)  , 

The  ejection  of  the  smoke  wet  studied  so  that  the  tip  of  • 
the  blade  cut  through  the  awoke  at  different  aslmuth 
angles  y  We  called  for  84  values  of  V  . 


Duo  to  the  fact  that  the  balance  ars  hindered  passage  of 

the  gmoke  for  the  obtainment  of-;0<T<  360*,  the  dlreo* 
tioo  of  rotation  was  changed  every  18  positions  and  noted 

ss  follows  :  •  » 

•  » 

*  * 

+ 'rotation  :  lBO'CT^  JdO* 


-  rotation  :  0<^<l80t  (see  photographs  4  and  8) 

Internal  Smoke  Emission  (Photographs  11  through  16) 

• 

This  form  of  emission  enabled  several  ^  pros  it  ions’  to  be 
grouped  on  a  single  photograph.  t 

e 

_ : _ l 
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there 

It  may  be  noted  (anc/will  be  further  reference  to  this 
when  discussing  the  data  processing  method)  that  any  one 
photograph  made  during,  say  a  +  rotation,  does  not  give 
all  the*  azimuths.  Indeed  it  can  be  seen  from  Figure  54 
that  a  number  of  blind  angles  arise  from  the  position  of 
.  V>e  cameras,  examples  in  +  rotation  being  340*< k>060*  and 
160*4  M'  4  180*.  The  cameras  were  therefore  moved  and 
■  further  photographs  made  to  supplement  the  results.  An 
*'  •  example  of  this  on  an  induced  velocity  curve  is  provided 

.*  ifl  Spftion  1.5  (cameras  at  C2  and  C5). 


the*  lhstantaneoue  position  pf  the  blade  was  also 
varied’.  In  general  ,  photographs  were  taken  of  the  oenter* 
plane,  advancing- blade  and  retreating-blade  positions. 

•  e 


Smoke,  Emission  Through  the  Hub  (Photographs  17  through  20) 

e 

a  •  * 

A  oentral  cooke  emlasAon  was  also  made  for  two  instantane¬ 
ous  blade  positions. 

•  • 

Zn  ftll  eases  and  for  each  form  of  smoke  emission  conside¬ 
red  above,  the  photographs  were  taken  from  vertical  and 
horltontal  viewpoints  respectively.  There  will  be  further 
oooaelon  to  discuen  the  difference  between  these  photo- 
gra phi  when' the  data  processing  method  is  dealt  with. 


in  eaoh  oate  a  recording  was  made  of  the  lift,  the  flappinb 

and  the  di*ag.  f 


102  ;-  Group  C  (Figures  |$  and  %6) 


In  group  0  only  films  •nd  recordings  were  made.  I 

The  teste  were  called  for  on  blade  No. 4  Tor  the  reasons  state 
previously  j  however,  this  blade  was  corroded  by  the  smoke 
.  froa  the  preceding  teste  and  was  therefore  replaced  by  the 
lighter  blada  flo;2  (0  .  0.093  a,  or  0.312  ft). 

'*  »  .  • 

The  Conditions  Of  deteralnatloa  of  t\ie  different  parameters 
were  "lore  rigid  than  preoedjcgly,  because  the. variations  of 
the  Input  variables  were  obtained  by  aeana  of  case,  and  thus 

effeoted  with  a  high  precision. 

♦  .  S''...  • 

The  input  variables  were,  the  sane  as  for  groups  A  and  B, i 

plus  the  tip  speed  U,  the  wind  velocity  ?©,  Two  extra  accele 

ration  tests  were  also  performed. 

•  • 

•Table  1.6  lists  the  test  program.  • 

.  and  only 

Smoke  emissions  were  of  the  internal  variety  on!y7one  case  of 
rotation  was  called  for.  However,  both  the  vertioal  and 
plan  views  were  retained.  I 


Each  test  was  accompanied  by  flapping  and  lift  recordings. 

/. 
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1,4  -  Data  Processing  Method 

The  processing  of  these  visualizations^  having  for  its 
object  to  determine  the  velocity  induced  over  the  diso  at 
different  azimuth  angles  V  ,  use  was  made  of  the  photographs 
of  groups  A  and  B,  in  a  different  manner  according  to  the  fori 
of  emission. 

1.41  -  Principle  Underlying  the  Method 

The  principle  however  is  identical  in  both  cases. 

Let  us  consider  two  vortices  (Figure  55&)  formed  at  the 
same  place  at  a  given  ^  by  each  blade  in  turn  one  at  an 
instant  Z  ,  the  other  at  T  +  Z  . 

2 

Let  these  vortices  be  designated  by  A  and  B  respectively. 

Relative  to  the  plane  of  the  rotor,  vortex  B  will  have 
been  carried  along  by  a  horizontal  velocity  V  and  by  a 
vertical  velocity  Vir  .  will  be  neither  the  induced 

velocity  at  B  nor  the  induced  velocity  in  the  swept  disc 
plane,  but  a  mean  induced  velocity  beneath  the  disc.  Vor¬ 
tex  A  will  have  sustained  a  same  entrainment,  but  during 
a  very  short  time  only;  furthermore  it  will  have  sustai¬ 
ned  only  a  mean  vertical  effect  Vix  close  to  the  disc  and 
of  lesser  magnitude  than  Viz .  1 

If  we  now  consider  A  and  B  relative  to  each  other,  we  may 
postulate  that  the  distance  AB  is  projected  horizontally 
and  vertically  to  give 

AA'  =  Vt 

AA"  =  Viz  t 

respectively,  where  t  is  the  time  elapsing  between  the  gem  - 
ration  of  A  and  B.  Since  A  and  B  are  formed  by  each  of 
the  blades  and  result  from  the  same  ^  ,  then  t  =  T  . 

2 

The  horizontal  velocity  V  will  represent  the  sum  of  thewlnd 
velocity  V0  and  of  a  velocity  Vix  (this  was  determined 
at  each  processing  performed). 

Once  this  principle  is  established,  only  the  modes  of 
applying  it  will  differ  according  as  to  whether  the  emis¬ 
sions  were  internal’  or  external.  These  differences  will 
in  fact  arise  especially  from  determination  of  the  angled, 

1.42  -  External  •Emission 

• 

Figure  55  groups  together  the  determination  of  (Figure 
55b,  horizontal  photograph)  and  the  determination  of  V< 
(Figure  55a,  vertical  photograph),  as  discussed  previuOis- 

ly. 
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Vortex  A  was  formed  Z  seconds  before  by  the  blade  A, 
and  was  formed  when  that  blade  out  through  the  smoke  fillei  . 
Let  A'  be  the  position  of  the  blade  at  the  instant  of  ge¬ 
neration  of  the  vortex  and  Y  its  azimuth. 

© 

1.43  -  Internal  Emission 

Contrarily  to  the  previous  case,  we  shall *choose  here  the 
angles  Y  in  order  to  locate  the  vortices  upon  their 
locus  (smoke  fillet).  Processing  of  these  photographs 
is  somewhat  tricky  ,  as  they  overlap  to  some  extent  there¬ 
by  entailing  numerous  correction  factors  regarding  the 
,  respective  positions  of  the  cameras  (Figure  56). 

Photograph  56b 

We  determine  the  angles  ((1)  (2)  (3)  )  along  the  trace 
of  the  swept  circle.  Neglecting  the  effect  of  Vy  ,  we  will 
therefore  have,  on  each  vortex  fillet,  the  position  of  the 
vortices  generated  at  (1)  (2)  (3),  say  (4)  (5)  (6)  for 
blade  B  and  (4  * )  (5f)  (6')  for  blade  A. 

Photograph  56a 

Viz  and  Vix  will  be  determined  from  this  photograph,  on 
which  it  will  suffice  to  locate  the  points  (4)  (5)  (6) 
and  (41)  (5*)  (61).  Figure  56  shows  the  construction  used 
for  the  purpose,  the  processing  taking  place  in  accordance 
with  thd  principle  indicated  previously  . 

In  both  processings  ,  and  as  regards  the  corrections  to  be 
made,  account  was  taken  of  the  respective  positions  of  the 
photographed  object  and  the  camera,  the  latter  having  in 
all  cases  been  focused  on  the  hub. 

Figure  5^  gives  a  schematic  diagram  of  the  relative  loca¬ 
tions  of  the  cameras. 

1.5  -  Test  Results 

From  a  study  of  these  tests  it  was  possible  to  deduce  the 
following  f ur.damer.tal  results  : 

-  the  evolution  of  the  mean  induced  velocity  beneath  the  disc, 
in  terms  of  the  flight  parameters  (photographs), 

-  the  experimental  lift,  drag  and  flapping  values  (recordings)  , 

-  the  assigning  of  circles  as  substitutes  for  the  helix  vor¬ 
tices. 

1.51  -  Evolution  of  the  Mean  Induced  Velocity 

Figures  57  through  60  give  the  rough  experimental  results 

/. 


for  two  selected  test  cases  computed  from  the  various 
types  of  visualization  defined  previously.. 

These  results  concern  the  Vc*.  and  Vix  values  respectively. 

The  following  flight  parameter  assumptions  were  made  for 
these  selected  test  cases  t 

ju  =  0.05  e^-  t*  ez-  o 

and  ® 

(ac  ©  ® 

P  =  o.io  ec=  i2*  <*«  3 ®  o  %■»  -  3 

Figures  60  through  64  represent  these  same  tests,  except 
that  the  curves  have  been  made  somewhat  smoother. 

The  evolution  of  Vtx  has  not  been  given  for  all.the  vlsuali  • 
zation  cases,  for  it  was  found  that  for  p  }  0.10, 

V1*  mean  = 

In  all  the  test  oases  we  find  evolutions,  in  the  case  of 
Vtj  ,  similar  to  that  of  Figures  61  through  63,  with,  in 
particular,  a  drop  in  Vi*  for  the  retreating  blade.  Figu¬ 
res  65  through  97  give  the  processed  results  for  various 
flight  parameters. 

An  8-point  Fourier  analysis  on  the  "typical"  Vt*.  curve  show  d 
that  the  third  harmonic  is  sufficient  to  provide  a  good 
approximation  (Figures  98,  99)  and  enabled  a  Vi'tm€4Jo  be 
calculated.  This  Vi%m  curve  is  drawn  for  different  flight 
conditions  (Figures  100  through  106),  as  is  also  the  non- 
dimensional  ratio  XiSs.  (Figures  107  through  113)  in  terms 
of  p  =  with  *  j/TST  ,  where  A«TrR*4ndp«^- 

and  where  Fjg  is  the  lift  given  by  the  recordings. 

The  following  equations  were  found  for  Vi_  * 

/*  -  0.05,  12*  3",  %  -  7\  V  0. 

Vcj,=  6.95  +  0.47  Sin  f  -  0.461  sin  2  f  -  O.557  sin  3^ 

+  0.745  cos  f  +  1.21  cos  2  Y  +  0.055  cos  3<K  (in  m/i  eo 

JJ  =  0.10,  ec=  12®,  =  3°,  0,  3* 

=  6.15  +  1.74  cos  V  +  1.37  cos  2  4'-  0.593  cos  3  ^ 

-  0.5  sin  2  V  -  1,2  sin  3  ^  (in  m/sec) 

In  these  two  test  cases,  the  VizB  values  are  6.95  m/s 'and 
6.15  m/8  respectively. 

for 

These  values  of  are  useful,  particularly  situating  the 
vortex  circles. 
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I.52  -  Experimental  Lift.  Drag  and  Flapping  Values 

&  9  9 

*  The  conventional  processing  performed  on  the  recordings  in 
each  test  case  enabled  the  lift,  drag  and  flapping  to  be 
defined  from* their  fundamental,  first  and  second  harmonics. 

In  Table  1.8  the  values  of  FN  ,  ?nx  ,  T0  ,  T,  ,  T{  and  a0, 
ai,  &2  are  given.  Also  added  Is  a  Table  1,9  giving  the 
rotor  lift*  in  ron-dimensional  fw  with 

1.55  -  Situating  the  vortex  Circles 

$ 

The  purpose  here  was  to  replace  the  vorte*  lines  (1)  (2) 

(3)  (Figure  114)  by  mean  circles  parallel  to  the  plane  of 
the  swept  disc. 


It  was  therefore  necessary  to  determine  the  distance  bet¬ 
ween  the  circles,  also  their  centers  and  their  radii. 


1.5)1  -  Distance  between  circles 


In  the  case  of  the  first  fillet  (1)  which  is  distant, 
in  time  ,  by  T  from  the  swept  disc,  a  good  approxi¬ 
mation  can  be  £  .  obtained  by  regarding  it  as  lying 
within  the  swept  disc. 


We  have  already  seen  that  the  fillets,  which  are'produ-| 
ced  in  alternation  by  each  blade,  are  separated  from 
one  another  by  a*  timfe  T  .  This  is  equally  valid  for 
the  circles,  s’ince  2  (1)  and  (2)  are  distant  by  T 
from  each  other,  then  if  (l)  is  located  upon  the  2 
swept  disc,  the  distance  between  (2)  and  the  swept  disk 
is,  in  mean  value.,  of. 3  T  .  Eut,  at  the  tip  of  the  blade 
this  distance  is  equal-’**’  to  T.  The  distance,  between  ijhe 
others  fillets, (2)  and  (3).. 2  ,  must  be  taken  as  T  . 

1.532  -  Circle  Centers  ^ 


All  the  circle  centers  are  remote  from  the  disc  plane; 
we  determined  this  experimentally  from  photographs. 


Ihesecenters  (Figure  114)  are  aligned  upon  a  straight 
line  ABC  suoh  that  0  A  «  V0  T  ,  where  V0  is  the 
relative  wind  velocity.  The  5  straight  line  ABC 
makes  an, angle  6  with  the  disc  diameter  lying  in  the 
center  plane  (  V  •  0,  V'’*  180*).'  We  first  determined 
0  experimentally  (Tables  I. 10,  I. 11);  we  then  oompared 
it  to  angles  CT  such  that  tanc/'»  Vtem^efc.From  tables 
1,10  and  I. 11  Jt  may  be  noted  that  J  .  This  result 

is  logical  enough,  for  the  straight  line  A  B  C  is  loca¬ 
ted  at  the  center  of  the  helix  and  is  not  subjected  to 
disturbance  by  vortices  whose  locus  is  the  helix,  or, 
if  it  is,  the  disturbances  caused  are  symmetrical  and 
cancel  one  another.  Hence  , 
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In  order  to  be  able  to  place  circles  in  the  hori¬ 
zontal  plane  we  have  defined  the  distance  OB'  (Figure 
114b),  where  B"  is  the  projection  of  the  center  of 
the  first  circle  on  the  horizontal  plane  passing 
through  the  center  of  the  swept  disc. 

0  B'  =  Vo  T  +  V(  J  +  \/i  g 


where  is  the  distance  between  two  ciroles 

tan  (3  -  Vo 

V0 

0  B'  =  V0  T  cos  +VoC0SK(nT) 

1 . 533  -  Circle  Radii 

It  will  be  remembered  that  the  circles  are  parallel  to 
the  .swept  disc,  so  that  they  form  an  angle  o<  with  the 
relative  wind.  t  *  * 

t  * 

It  now  remains  to  determine  their  radii. 

To  this  end  *we  shall  first  determine  experimentally  the 
directions  and  eJ*.  defining  a  bounded  zone  containing  th 
circles  (Figure  114).  In  the  zone  >  180*,  the  direc¬ 
tion  o4  does  not  bear  upon  the  blade  tip  ,  but  is  at 
a  distanoe  O'A  ■  V0  T  . 

y  * 

From  the  photographs  we  found  different  values  for  €*  tf, 
o\  for  0  <.  0.05.  On  the  other  hand,  for  u  ^O.IO,* 

we  have.  ' 
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to  w 


The  values  f  ^  $  are  capable  of  explanation.  Let 
us  call  Vi.*m  the  mean  horizontal  induced  velocity  (Fi¬ 
gure  114b).  The  vortex  at  E  will  be  influenced  by  the 
Vix  of  the  vortex  E'  and  will  be  accelerated  with  respec 
to  the  relative  wind.  On  the  other  hand,  the  vortex  F' 
will  be  retarded  by  the  influence  of  F’. 

We  can  therefore  calculate  ^  and  <5^  from  the  following* 

tan  & .  tan  S,  - 

V<l*W*in  vo-Kxm 


As  already  stated,  is  defined  from  the  data  prooes 

sing  in  each  test  case.  Also,  when  V0  beoomes  relative¬ 
ly  large  (p>  0.10),  we  obtain  Vi  ^0,  which  explains 
the  equality  between  c?  »  .  <SZ  . 

The  radius  of  the  circles  is  given  by 


Ass  R  __  14-4  *  Vcz  S  tecfj  ^ 


-A-33- 


<r>  % 

°  *  • 

th» 

where  Viz  £  is  distance  between  two  circles 

Ro  is  the  radius  of  the  swept  disc  and 

tan  and  tan  have  been  defined  previously. 
Let  us  take,  for  0  <  ym  <  0.05, 

R_  R  V'T+i  Vix*Z 

0  z 

and  for  0. 10  <  jU  <0.20 

R  =  R0  -  Va  t 
2  IT 


In  Tables  I. 10  and  I. 11  we  summarize  ,  for  different 
parameter  values,  the  values  of  V0,  Vi,.  ,  Vit_M  ,  of  S 

calculated  ,  of  u  measured  and  of  Vixra  (  for  group  B) . 
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BASICS  FOR  THE  ESTABLISHMENT  OF  AN  ANALYTICAL 
METHOD  OP  EVALUATING  AERODYNAMIC  LOADS 


2.1  -  Causes  of->a  discrepancy  between  the  predicted 

d,  lift  and  the  actual  lift  of  a  blade  element 

* 

2.11  -  The  local  lift  dL  of  a  blade  element  c  dr  (see  Figure 
.  115)  Is  usually  calculated  from  the  following  formula: 


(II. 1. 


The  values  of  the  flapping  coefficients: 


(II. 1.2)  /*>s40  (t)-«3^(t)  cos  f-daCt)Co5  IV-A3  <\)  CoS  3^ 

_  bH(t)3in^-  bz  (k)  sin  ZV -  b3  (t)  sin  3  Y 


(II. 1.3) 


are  calculated  (hinged  blade)  from  the  equation  : 


f  *L 
7  ^ 


dJz.  r  cii 


(II. 1.4) 


The  induced  velocity  V*  is  in  general  defined  by  overall 
formulae.  The  simplest  of  these  is  : 

y.  _  _ K 

c  zf  a/v*^v?. 

Discrepancies  can  arise  for  a  variety  of  reasons. 


2.12  -  Simplification  of  the  formulae 


In  actual  fact,  the  local  lfft  should  be  expressed  by  : 


(II. 1.5) 


There  is  an  error  in  the  velocity  t 


4-Vsin  Y)  +  (V^  +  VsincL+r  ^/JcoS^KJ 


(the  second  square  term  is  often  neglected). 
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«  There  is  an  error  in  the  evaluation  of  the  angle  of 
!  Incidence  : 


(11. 1.6)  . 

(11. 1.7) 


t  c  0-4* 


0  «  0C  +  9^  cos  V  ♦  &z  sin  ^ 


(II. 1.8)  _  Vi+Vs\n*+  rff  4-VccS^t  fc.CoS^ 

- <or  +  Vjte  V 

The  tangent  arc  and  the  angle  are  taken  to  be  one  and 
the  same. 

This  gives  rise  to  an  error  of  non-linearity;  if  the 
angle  of  incidence  i  becomes  too  large  it  is  no  longer  poss 
to  replace  Cl  by  (  9-  f) . 


If  one  confines  oneself  to  low  advance  ratios 

w  =  -4-  4  o.z 

>  cu  R  ^ 


and  if  one  avoids  over-accentuated  stalls,  one  may  retain 
the  formulae  (II. 1.1)  and  (II. 1.5). 

2.15  -  Blade  rigidity 

Blades  are  usually  flexible  in  bending  and  sometimes  also 
not  very  rigid  in  torsion; 

For  a  flexible  blade  fi  no  longer  has  any  physical  signi¬ 
ficance.  If  y  (r,t)  be  the  distortion  of  the  blade  (see 
Figure  116)  the  expression  for  Ji  must  be  replaced  by  the 
following  expressions  : 


e> 

by 

Jr  * 

a  3 

by 

3  u 

dt 

3  b  d 

d*fr 

by 

a‘.M 

dbz 

at4  e 

The  blades  in  the  five  sets  built  were  made  very  rigid 
in  order  to  avoid  the  y  distortions  at  low  frequencies. 

A  parasite  rotation  of  the  blade  can  also  take  place 
if  the  blade  distorts  in  torsion  or  if  the  pitch  controls 
.  are  "mushy". 

The  blades  built  for  the  te<ts  were  torsionally  rigid. 


* 


» 


t 
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The  CG's  of  the  various  blade  sections  were  located  upon 
the  feathering  axis.  • 


On  the  other  hand,  the  controls  were  light  enough  to  per¬ 
mit  remote  control  and  to  introduce  unit-step  harmonic 
f  orces . 


It  is  considered  that  the  parasite  cf©  effects  from  the  • 
inserted  blade-roots  were  small. 

Consequently* the  expressions  (II. 1.1)  and  (II. 1.2)  remain 
valid. 

« 


2.14 


Induced  velocity  and  blade  flapping 

A  priori,  the  only  doubtful  terms  in  equation 

are  :  , 

*  * 


V 


o  > 


Aft  .  „ 

r  -  and  ft 

At 


♦ 

(II. 1.1) 


The  induced  velocity  Vj  is  generated  by  the  free  vortices 
escaping  from  the  t  rotor  blades  (b  =  2). 

As  long  as  the  free  vortex  is  remote  from  the  blade  (at 
a  distance  of  about  3  to  4  chords)  it  produces,  on  the 
straight  section^  of  the  airfoil,  a  constant  velocity 
distribution  that  generates  an  ordinary  type  of  circula¬ 
tion  (see  below):  hence  the  blade  can  be  reduced  to  a  line 
and  the  velocity  Vi  sought  at  that'point. 

When  the  vortex  is  close  to  the  blade  (at  a  distance'of 
less  than  3  chords) ,  the  blade  can  no  longer  be  reduced  to 
a  straight  line  because  the  distribution  of  the  velocity 
induced  by  the  vortex  on  the  airfoil  varies  too  much, 
thereby  causing  inertia  forces  to  be  set  up  in  addition  to 
the  circulations. 


The  variation  through  time  of  the  'Circulation  about  the 
blade  causes  free  vortices  to  escape  from  the  trailing  edg<  . 

These  two  latter  free  vortex  effects  produce  dynamic  (non¬ 
stationary)  regimes  that  must  be  taken  into  account.  Thes( 
will  be  examined  again  later.  As  a  result,  the  lift  dL 
may  sustain  amplitude  modifications  as  well  as  an  out-  S  ^ 
phasing. 

The  static  lift  dLstat>  given  by  formula  (II. 1.1)  conse¬ 
quently  becomes  a  dynamic  (non-stationary)  lift  d£  dyn 

vdr  dyn.  *  '3?  '  st*t. 

where  designates  an  operator,  i.e.*  a  mathematical 

_ u _ 
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(V 

transformation  of  the  static  lift,which  takes  into  accoint 
the  amplitude  and  phase  modifications. 

The  blade  flapping  value  is  consequently  modified: 


(II. 1.5)'  + 

2.15  -  In  what  follows, 
and  the  operator 


it  is  proposed  to  attempt  to  show.how  Vi 
can  be  calculated.  ' 


2.2  -  Description  of  vortex  model 

2.21  -  General 

» 

Before  proposing  a  calculation  method,  a  vortex  model  will 
be  described  and  the  orders  of  magnitude  of  the  various  phe¬ 
nomena  estimated  §0  as  to  permit  simplification  of  the  ma¬ 
thematical  equations. 

• 

In  the  case  of  an  aircraft  fixed  wing(in  a  compressible 
medium),  the  whole  set  of  phenomena  is  known.  On  the  other  hand 
whereas  the  phenomena  are  similar  in  the  case  of  a  helicop¬ 
ter  rotor,  expressing  them  mathematically  is  a  complex  busi¬ 
ness;  an  attempt  will  therefore  be  made  to  find  approxima¬ 
tions  for  these  mathematical  expressions,  by  introducing 
the  phenomena  in  rough  form  only. 

The  vortex  visualizations  obtained  by  means  of  smoke  emis¬ 
sions  proved  of  great  help  in  achieving  these  approximations , 
and  especially  as  a  means  of  insuring  that  all  effects  had 
been  accounted  .for. 

2.22  -  Helix  vortices  and  radial  vortices 

The  free  vortices  are  generated  on  the  blade  and  are  of 
two  kinds. 

Any  change  in  the  circulation  P (  r,t  )around  the  blade  : 

(II. 2.1)  dL  e  p  V,w  r  dr 


gives  rise  to  a  streaming  of  free  vortices. 

2.221  -  First,  the  circulation  H  can  vary  along  the  span  t 

.  dr 

These  free  vortices  locate  themselves  along  the  flow 
streamlines  if  they  are  referred  to  axes  linked  to  the 
blades 


(  V  .  rotation  V  =  0  ) 
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These  .pill  be  germed  helix  vortiees  (Figure  117). 

Let  us  now  calculate  a  point  of  such  an  helix  vortex 
(Figure4ll8) . 

Let  us  consider  the  blade  at  the  instant  t  -  ti  and  a  sys¬ 
tem  of  coordinates  0  x  y  z  linked  to  the  blade,  and  let  M 
be  a  point  on  the  blade  from  which  is  detached  an  helix 
vortex  of  intensity 

The  coordinates  of  the  point  M  are  t 


I  r 

(II. 2. 2)  <  lK(tO 

l  1*  tM 


{acM  s  -rtos 
u  M  ‘  s  r  s\f\ 


Let  us  now  find  the  location  of  the  helix  vortex  element 
formed  Z  seconds  before,  i.e.  at  the  instant  t> . 


Z  seconds  before,  the  center  0  of  the  rotor  was  at  0^2  • 


*  "vf  Z 


(II. 2. 3) 


Z  seconds  before,  the  point  M  was  at  M2,  suoh  that: 
x  u  s  -  /*  cos  tt>  ( b.,  -Z)  -  Vc eS  ok  Z 
^  s  /•  sin  wu  -Z) 

Ct1-  Z)  -  Vsino*S 

At  that  instant  (  fcA-  ®)  a  free  vortex  element  became  deta¬ 
ched  at  Mg  . 


During  the  time  taken  for  the  blade  O2  M2  to  reaoh  0  M 
(viz.  two  seconds)  this  vortex  element  was  entrained  from 
the  point  M2  to  the  point  H  by  the  velocities  induced  by  al 
the  free  and  connected  vortices  t 


* 

vc* 

<■  .2 

(11.2,4) 

c  *■ 3 « ; 

&*l) 

*5 

1  x'  3 

> ,t; 

The  coordinates  of  H  are  therefore  t  % 

c-rcosu;  (  t^-S)  -\/cos  Z  +  £  VL^AZ 

u 

hi 


r 

■  r  jinw  C  -Z)  e*  a  d  Z 

Jo  *  z 

■  ~  r  (3  -  Z )  +■ 


(II. 2. 5) 
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The  slope  of  the  vortex  line  at*H  is  consequently: 

# 

**  m  -  r  (ju  Sin  w  (t  4  -  S  j  -  Vco*  d.  +  V/cx 

d  Z 

AOjl  »  -  r  u>  cos  co  ( -  $)  +  N/ca  . 
cLZ  * 

d^H  -  r  £H  (>i  -S)  +  yjirt*  +v/i,» 
d£  dfc 

Let  us  consider  the  difference  with  an  aircraft  fixed 
wing. 

In  the  case  of  fixed  wings  it  is  customary  to  cancel  out 
Vi*  ’»  v<-4  »  V;.*,  in  the  wake,  since  their  values  are  small 
compared  with  the  flying  speed  V0. 

From  this  it  may  be  asserted  that  the  free  vortices  remain 
at  the  place  (referred  to  space  coordinates)  where  they 
were  formed  (f^)* 

This  is  not  always  true  in  the  case  of  a  helicopter. 

The  helix  vortex  is  tangential  to  the  streamline  at  that 
point. 

Determination  of  ,  VLu  and  of  is  a  difficult  problem 

The  visualization  tests  made  have  nevertheless  enabled  vali 
approximations  to  be  made. 

2.222  -  Second,  the  circulation  P  may  vary  at  a  point  r  on  the 
blade  in  terras  of  time,  as  follows  : 

bV 

at 

•,/hen  that  is  the  case,  a  vortex  will  stream  from  the  trai¬ 
ling  edge.  Such  a  vortex  will  be  termed  .a  radial  vortex 
(Figure  119)  and  its  intensity  will  be  equal  to  : 

_  i£_ 

dfc 

since  it  is  opposite  to  the  intensity  which  generates  it 
(conservation  of  the  circulation). 

• 

This  radial  vortex  is  orthogonal  to  the  streamlines  defi¬ 
ned  above 

V  TT 

(  V  .  rot.  V  =  maximum  )  —  -j- 


The  coordinates  at  the  instant  t  ■  ti  of  a  radial  vortex 
•  / 
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(II. 2. 7) 


element  formed  %  seconds  before,  l.e.  at  the  Instant 
ti  -  Z  ,  are  j 

rz 

xr  s-r  cos  ou  (t*  -Z)-VootoLt  Vix'dZ 

g  * 

.7)  *  r  s,rx  00  VCy  dZ 

JT  s  -rj*  (t4  -Z)  +  VSiOotZ+J*  WjfllS 

•  The  tangent  to  the  radial  vortex  line  is  defined  by  : 


^*7 


-  -COS  CD  + 


Vl  xdZ 


(II. 2. S)  .  sin  w  +37-/  ^dZ 

*^f0  VildZ 

The  orthogonality  of  the  radial  vortices  and  of  the.  helix 
vortices  can  be  expressed  by  : 

iir .  it*a  +  iix.iiS.  +  4^-  4^  -° 

IT  X  +  17  «  1  Tt 

2.23  Conservation  of  circulation 

2.231  -  The  complete  system  of  free  and  connected  vortices  consti¬ 
tutes  a  system  of  closed  curves  of  elemental  intensity 
(Figure  121): 

a^r 

Each  elemental  closed  curve  consists  of  connected  vortex 
seotion  (otft),.an  helix,  vortex  section  (  (1  ?) ,  a  radial 
vortex  section  (  )  and  a  further  helix  vortex  section 

(  )  • 

•  The  free  vortex, viz. (<£  excluded)may  assume  a  wide  variet; 
of  configurations  and  change  its  location  periodically  inn 
ponse  to  the  velocity  field  (depending  on  whether  operation 
takes  place  as  a  helicopter,  a  gyrodyne,  or  in  autorota- 
ting  descent);  the  total  circulation  will  consequently 
by  zero.  It  is  proposed  to  consider  only  helicopter 
,  configurations  here. 
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Por  long-established  hovering  ($-**•)>  the  radial  vortex  tfo 
is  thrown  back  to  downstream  infinity  ( &  2»  )  and  th< 
helix  vortices  and  cJot  are  of  infinite  length.  1 

2.232  -  The  connected  circulation  produces  the  local  lifti 

P  P  (  r,b) 

Since  the  lift  cancels  out ‘at  the  two  ends  of  a  blade  (at 
the  root  (V0 )  and  at  the  tip  (R)  (see  Figure  120)),  we 
have  j 

T(  r0 , b)  *0 

r  c  ,  t  >  *0  .. 

* 

The  conservation  of  circulation  gives  : 

♦ 

n 

(II. 2. 9)  /  iJL  (  r,fc  )  dr  «  r  (K)  cr0)  .0 

Jr  5  p 


Whence  : 


/  r* 

P(r)  •  /  (r,t)  ir  = -/ 

J  dr  Jr  dr 

'4 


The  circulation  at  a  point  in  the  helix  vortex  is  \ 


AL  (r^rzJhL 

dP  1  1  /  OrdZ 


dz 


The  circulation  at  a  point  in  the  radial  vortex  is  : 


ar 
a  t 


r  1 


£SL  aLf" 

e  ara  & 


2.233  -  Experimental  orders  of  magnitude 

In  order  to  estimate  the  values  of  : 


/>,<n,  $£(r,1)rt  (r,f) 


these  will  be  calculated  from  the  experimental  lift  dis- 

/. 


(II. 


V 


4 

tribution  values  £i  Min)  , based  on  the  tests  carried  out 
by  John  Meyer  Jr.  and  Gaetano  Falabella  (see  refe¬ 

rence  .4.). 

•  Tv/o  cases  will  be  considered  :  hovering  flight  (  u  ■  0) 

and  forward  flight  (^4  ^  0  ).  The  rotor  characteristics  ars 
given  by  Table  2.13  and  the  circulation  P  is  obtained 
from  : 


P  A  dL 

“  F+  fusing  dr 

The  points  in  the  neighborhood  of  r  =  R  and  r  =  r0  are 
somewhat  uncertain. 

The  variations  of  P  in  terms  of  r  are  given  by  the  graphs 
in  Figures  122,  123,  124. 

The  maximum  value  of  P  is  situated  in  the  vicinity  of 
r  =  0.8,  for  =  0.2  and  r  '=  0.9  for  p  =  o. 

The  values 


P  (  r  *  0.5,  ^  )  and  P (  r  =  0.8 ,  Y  ) 


have  been  broken  down  into  a  Fourier  series  (Runge's  12- 
point  method)  for  forward  flight  ( >*.  »  0.20). 

p 

The  expressions  for  —  are  given  by  Table  2.12). 


.10) 


[  rn  COS  n  V  + 


sin  nV 


As  a  first  approximation,  we  may  write  : 


This  result  bears  a  relation  to  that  found  by  Meijer  Drees 
(see  reference  Jf, ) , 


l'f  it  is  assumed  that  the  circulation  f*  (r) 
the  span, •  so  that  t 


.  0 


dr 
TF 

then  the  blade  flapping  equation  j 


is  constant  over 


(II. 2. 11) 


(aA  *  V  c uzli  +•  1.1*  « 
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developed  in  cosine  and  sine  form  gives  the  following  two 
simple  relations  : 


If  the  hinge  is  a  central  one  {a  =  0)  then  we  have  : 

n/ 


n> 


and 

£-j  _ 3  u 

I Y  '  f  h 

Using  the  local  values  for  T  given  by  Table  2.12,  it  will 
be  seen  that  the  second  relation'  : 

■A*  -  -  0.15 


mr 

is  true  neither  in  sign  nor  in  magnitude,  for  : 

P4  f  +  1.47  for  r  *  0.5 

^  +  1.52  for  ?  =  0.8. 


On  the  other  hand,  the  first  relation  : 

r»*  Pi 

•P#  p*  P»P 

gives  0.415  and  0.245  for  with  -—4-  negligible. 

From  this  it  emerges  that  the  hypothesis  of  constant  cir¬ 
culation  is  untrue  for  cos  ^  and  sine  2  V  frequencies  and 
that  it  would  provide  an  approximation  for  the  0  and  Sia  V 
frequencies,  unless  the  .blade  flapping  equation  (II. 2. 11)  is 
incomplete , 

dP 

was  calculated  by  deriving  the  curve  T  (r)  graphically 
This  graphical  derivation  was  performed  as  follows  : 

a)  A  protractor  was  used  to  measure  the  slope  of  the 
tangent  in  degrees.*  . 

b)  A  curve  was  drawn  to  represent  the  values  of  the 
angles  found;  this  curve  was  then  'smoothed  out' 

/  * 
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Table  2.12 

Harmonic  Analysis  of  the  Circulation 
around  a  Blade  at  f  ■  0.5  and  r  *  0.8 

for  jA -■  0.20  (Sec  Reference  ...) 
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2.234’ 


and  checked  against  the  measurements  obtained 
under  a)  above. 

c)  The  angles  were  converted  into  trigonometry  tan¬ 
gents,  thus  enabling  the  derivatives  to  be  obtai¬ 
ned. 

The  variations  jpr  (  r,Y  )  are  shown  by  the  graphs 
in  Figures  124,125  and  126. 

The  area  comprised  between  the  curve  and  the  x-axis 
must  be  zero  : 

r»  ft 

AIL  d t  «o 

a  r 

ar 

Although  the  value  *£7:  appears  to  tend  towards 
infinity  (00)  as  r— *.R  and  r— . *r0,  the  area  bounded 
by  the  curve  and  the  r  axis  between 


L 

'  *  ~o 


p  =  ''m 

(  r  5  rm  , k  )  -0 

and  R,  must  have  a  finite  value  ; 

r-*R 

dr 


4L  (r.Y)dn  *r 
d  r* 


max 


r*  r, 


m 


Determination  of  the  Induced  velocity 
as  a  complementary  term. 

If  it  is  assumed  that  the  local  circulation  is  expressed  byi 

d  L 


J-  =  p  (Ujr  +  V  sin  Y)  r 


r.  1  c  f«,r+  Vain  fifitl  f e  V^  +  l'  di  VcosYW/jind. 
1  U<JL  u>r+  V/sin'M' 


then  the  induced  velocity  Vi  will  be  equal  to  t 


S  U  |V  4 p.  srn  Yj  ©-  jr  ^2.  +  VcoS  Vj»nolj  - 


r 


(II-2.12) 
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In  the  case  of  the  Meyer-Palabella  tests  (see  referenoe  ... 
and  Table  2.13)  all  the  terms  on  the  right-hand  side  of  the 
equation  being  known,  it  was  possible  to  calculate  the 
induced  velocity  required  to  balance  the  equation.  Evi¬ 
dently,  all  errors  and  parasite  effects  are  included. 

The  Vi  curves  (  **  ,  V )  are  given  in  Figures  124,  127,  and  128 
for  hovering  flight  and  al30  for  forward  flight  (/*=  0.20). 

These  curves  were  compared  to  the  mean  induced  velocities* 


M3) 


for  0 


2f7rRVv/+Vi* 


for  fj  =  0.2 


These  calculations  call  for  the  following  observations  x 

1)  In  forward  flight  ( p  f  0),  the  induced  velocity  is  not 
Just  an  indirect  measurement  of  the  lift.  For  V  =  90°, 
due  account  was  taken  of  a  possible  error  in  G  of  0.5° 
(dotted  curve);  however,  it  does  not  account  for  the 

pr ssibility  of  the  negative  induced  velocity.  The  latter 
is  presumably  due  to  the  free  vortices  of  the  second 
blrde,  which  are  located  beneath  the  rotor  disc  (see 
visualization  photographs  ). 

2)  This  way  of  determining  the  induced  velocity,  from  the 
Meyer-Falabella  results,  i3  not  physically  correct. 

For  it  is  not  the  incidence  but  the  chord  which  cancels 
ouf  at  the  blade  tip,  thereby  giving  rise  to  a  large 
increase  in  the  circulation  derivative. 


/. 
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T&ble  2.13  m 

Data  Concerning  the  Meyer-Falabella  Rotor 


Diameter 

Airfoil 

Chord 

RPM 

Angular  velocity 
Central  hinge 

Hovering  flight  ^1=0  . 

Sc  = 

T 

Forward  flight  p  =  0.22 

9c- 

P  = 

CT  _ 

<r 

ot  = 


*2  R  -  5  ft 

NACA  0015  *  5.7 

d.  v 

c  =5  inches 
N  =  800 
CO  =85.8  rad/sec 

6  =  0 

8“ 

7.6  lbs 

8* 

0*5  -  4*6  cos  'K  -  2*  sin  V 
0.085 


5 
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2.24  -  Vortex  visualization 


-  The  accompanying  photographs  taken  at  the  Marseilles 
Wind-Tunnel  show: 

-  the  extreme  marginal  vortex  and 

-  the  vortex  at  the  blade  root. 

-  The  vortex  at  the  blade  tip  is  clearly  localized:  it  is 
also  more  marked  (sharper  photographs)*  for  blades  locate< 
between  azimuth  angles  »  180°  and  }60#  (retreating  bla< 
*des)  than  for  blades  located  at  0  <(fj<  180°  (advancing 
blades) (photographs  markedly  less  contrasted). 

core 

-  There  is  a  single  /6nly  in  the  cas§  of  a  steady  periodic 
regime.  According  to  a  film  which  was  taken, if  sudden 
collective  pitch  changes  are  made,  small  cores  may  appeal 
which,  after  a  certain  time  has  elapsed,  are  sucked  in 
by  the  main  core. 

-  The  vortex  at  the  root  is  much  less  intense  and  difficul' 
to  place  in  evidence.  The  spirals  correspond  to  a  high 
local  f*.  (  r  =  0.2) . 

2.25  -  Calculated  orders  of  magnitude 

2.251  -  Application  to  the  two-blade  rotor 

In  a  manner  somewhat  analogous  to  the  theory  of  fixed 
wings,  it  is  necessary,  in  order  to  calculate  the  lift 
at  a  point  M  on  a  blade  A,  to  know  the  velocities  indu¬ 
ced  there  by  the  free  vortices  and  the  connected  vor¬ 
tices  . 

It  is  therefore  necessary  to  kno\*  both  the  location 
and  the  intensity  of  the  vortices.  Visualizations  permit 
a  simplification  of  the  problem  of  situating  these 
vortices . 

Lpt  M  be  a  point  where  we  wish  to  know  the  velocity  in¬ 
duced  by  the  vortices  : 

=  -  r  coo  f 

(II. 2.14)  =  +  r  sin  $ 

= • "  r 

To  calculate  the  induced  velocities,  use  will  be  made 
of  the  Biot  and  Savart  formula  written  in  Cartesian  eooq- 
dinates  (Figure  129). 

This  point  M  will  be  on  blade  A  if 

f  -  % 

_ fin  -  _ A 
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It  will  be  on  blade  B  if 


^  a  ^  ♦IT- 

s  |i  (  )  where  T=period  of  rotation. 

2.252  -  The  coordinates  of  a  point  H  in  the  helix  vortex  are 

known  in  terms  of  the  auxiliary  variable  (see  Eq.II.2.5). 

Formula  II. 2.16  given  on  the  following  page  permits  calcu¬ 
lation  of  the  velocity  induced  at  a  point  M  by  the  helix 
vortex  (  Z  varies  between  0  and  00  ) . 

Since  it  was  necessary  to  make  approximations,  a  number 
of  numerical  applications  were  made  by  graphical  integra¬ 
tion. 

2 . 253  -  Helix  vortex  iJi  hovering  flight 

The  velocity  induced  at  a  {Joint  r  =  0.5  R  on  the  blade  A 
by  an  helix  vortex  of  intensity  P  streaming  from  the 
point  P^  =  0.8  R  with  : 

U  =  Oo R  =  100  m/sec 

R  =  O.75  m 


Vi  =  4.55  m/sec. 


is  given  by  : 


V-  =  — 
u  47T 


L  A* 


(r,  -  P,  r  toS  ujZ)  d  (ouZ*  ) 


(II. 2. 15) 


l  R  +f't-Zrci  +  ^ 

The  integration  is  performed  grapnically  for  the  follo¬ 
wing  three  cases  (Figure  130)  j 

1)  For  the  helix  ,Vi  =  4.55  m/s  0  ^  1 uZ  ^  2  n  ) 

2)  For  a  circle  (Vi  =•  0)  contained  in  the  plane  of  the 
disc  and  acting  as  a  substitue  for  the  helix. 


3)  For  a  circle  located  half  a  period  beneath  the  disc,  i.e, 


<A.)VuC)2=  (—)l  n-4 

'uR  '  lU)R/ 


7^Cf& rtcj Ls*  re  -  v^^oc/t-y  W7“^?  f*&srr7-  s-r  f  r*j  ^ 

&  Y  -4/V  |/i£?>^^£LV  AfJ4Avy$  WO/¥  ^  #  ^ 
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For  a  whole  period,  i.e,  for  one  complete  loop,  the 
following  values  are  respectively  obtained  for  A  (  1  ) 

(  meter) 


A  .  1  ® 

meter 

error  for  one  period 

1 

9.05 

0 

0  % 

2 

10.7 

+  .1.65 

18  % 

3 

9.7 

+  0.65 

7  % 

9 

For  a  quarter-period  ,  l.e.  for  the  first  quarter-loop  , 

we  have 


A.  1 

meter 

error  or  th 

e  first  1/4 

1 

4.5 

0 

0 

2 

4.5 

0 

’  0 

3 

4.0 

-  0.5 

-  12  % 

If  the  quarter  helix  be  replaced  by  a  straight  line  half- 
vortex  Issuing  from  the  same  point  (0.8  R)  and  having 
the  same  intensity,  we  have 

V,  _  JL  .  - i -  [cos  f 

4  rr  (o. 8-o.s;r  |_ 

f  — w.  o 

The  value  of  A'  then  becomes  : 

A'  =  4.45  m"1 
From  the  calculations  it  can  be  asserted  that  : 

a)  The  effect  of  a  quarter  helix  upon  its  own  blade  can  be 
replaced  either  by  a  quarter  circle  contained  in  the 
plajie  of  the  disc  or  by  a  straight  half-vortex 
(Figure  lb). 

b)  Tha  effect  of  a  complete  helix  loop  can  be  replaced  by 

a  circle  located  '  -.nealh  the  swept  diac,  one-half  period 
from  it  .  From  the  processing  made  of  the  visualiza¬ 
tions,  one  is  led  to  si  uating  the  circles  in  the  manne 
indicated  below,  in  order  to  take  into  account  the  de¬ 
formation  sustained  oy  the  vortices  :  / 
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-  When  calculating  the  effects  of  the  helix  vortices  on 
the  blade  A,  the  effect  of  the  latter's  marginal  helix 
vortex  can  be  replaced  by  that  of  a  circle  located  a 
quarter  period  beneath  the  disc. 

e 

-  To  study  the  effect  of  the  first  half  of  the  loop,  the 
half-loop  in  question  can  be  replaced  by  a  semi-circle 
situated  in  the  plane  of  the  disc.  This  gives  a  good 
approximation. 

-  To  study  the  effect  of  the  helix  vortex  of  blade  B  on 

blade  A,  the  effect  of  the  first  loop  can  be  replaced  by 
that  of  a  circle  located  3  T  beneath  the  swept  disc 
(Figure  131a) .  4 

The  effect  of  the  remaining  helix  loops  can  be  replaced  by 
circles  located  at  distances  measured  in  half-periods . . 

2.254  -  Velocity  induced  on  the  blade  A  by  the  first  loops  of  the 
marginal  helix  vortices  from  blade  A  and  blade  B. 


The  circulation  of  these  vortices  is  constant.  The  graphics 
integration  was  performed  for  forward  flight  (u=  0.20)  in 
the  following  two  cases  t  '  ' 


Case  (a)  1  Blade  A  is  in  the  position  *V  =  0  and  the  induced  v 
city  is  calculated  spanwise. 


bln 


Case  (b)  ;  Blade  A  occupies  different  azimuthal  positions, 
the  induced  velocity  being  calculated  for  a  fixed  point  P 
(r  =  0.75  R). 


Characteristics  are  as  foil ows : 


b  =  2  blades 

"R  =  0.7^6  meter 
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Case  (a) 

The  induced  velocity  at  the  point  */=  0  is  given  by  the 
table  below  : 

Table'  2.14 


r 

rs  TT 

A  V-A 

avc.,5 

ToUt  \Ji 

- 

"Vs 

m/s 

.m/s 

0 

0 

0.4 

1.0 

0.75 

0.625 

0.555 

1.18 

0.95 

•  2.-0 

5.0 

5.0 

It  may  be  noted  that  the  effects  arising  from  the  two  blades 
are  substantially  be  same  (  =  0). 

Case  (b) 

Velocity  induced  at  the  point  r  =  0.75  R  with  'f'*  assuming 
16  different  positions. 

The  result  is  given  in  the  table  on  the  following  page  and  is 
illustrated  graphically  in  Figure  1J2. 

It  may  be  noted  that  the  constant  -circulation  helix  give 
variable  induced  velocities. 

The  effect  produced  by  the  helix  of  blade  A  is  always 
sizeable . 


The  effect  of  the  helix  from  blade  B  is  particularly  marked 
in  the  second  rear  half  of  the  swept  circle. 


2.255  -  Practical  estimation  of 


Number  of  helix  periods 
for  the  following  case  : 


infinity 

necessary.  This  study  was  made 

(constant) 

0.C765 


b-  2 

Q  =  0.927  sq.m/sec 
R  =  0.746  ra 


fA=  0.2 

U>  =  111  rad/sec 


A 

o( 


I  t  «  I 
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A  graphical  calculation  was  made  of  the  velocity  induced 
at  the  point  r  =  0.75  and  f  i  0  by  the  helix  vortices 
up  to  the  second  loop;  beyond  this,  the  asymptotic  formula 
deduced  from  the  formula  (II.2.16)  given  previously  was 
adopted.  , 

A  ~  f  ut+A*  Cos*(vZ  d(<vZ) 


Q 


i-trH  (n*-+  A*-)  & 


C  tu  Z ) 


.00 


The  numerical  results  obtained  are  given  in  the  following 
table  2,l6. 


Table  2.16 


The  first  two  periods  give  the  induced  velocity  with  an 
approximation  of  4 

Le^  us  now  seek  the  location  of  the  downstream  infinity  for 

ft  =  0.2. 

The  first  period  gives  77.5  of  the  total  induced  velocity. 

The  first  two  periods  give  the  induced  velocity  at  4  JS  near 
of  the  total  induced  velocity. 

Therefore  the  second  loop  must  be  located  at  a  distance 
(Figure  151c)  : 

A0  A2  =  VQ  x  2  T 
where  T  is  the  period. 

We  have  : 


P  ,  o.2 


ouR 


T  = 


zn 

(V 
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whence : 

*Aq  A2  =  1.5  .  2  R 

The  practical  infinity  can  consequently  be  defined  as 
follows  : 

1)  Depending  on  the  value  of  the  advance  ratio  and  upon 
the  disc  loading  A  ,  n  helix  periods  will  be  taften 
such  that  : 


2)  The  practical  downstream  infinity  is  located  at  a 
distance  of  1.5  times  the  rotor  diameter. 

This  result  can  be  compared  either  to  half  a  straight 
line  vortex  (Figure  133)  or  to  a  solenoid. 

2.26  -  Conclusion 

These  basic  considerations  will  be  used  in  approximate 
methocb  of  calculation  exposed  in  the  following  chapter. 
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3  -  PROPOSED  METHOD  OF  EVALUATING  THE  AERODYNAMIC 
LIFT  OF  A  HELICOPTER  BLADE 

3.1  -  Hovering  • 

3.H  -  Setting  In  equation  form 

Taking  the  case  of  a  two-bladed  rotor,  the  local 
circulation  (r)  at  a  given  point  on  the  blade  A 
is  provided  by  the  expression 


(III. 1.17) 


H  =  -Icuf'c 

A(r)  Z  di 


i  L  wr  cur 


where  Vi®^  and  are  the  velocities  induced  on  a 
station  point  or  blade  A  by  the  free  vortices  issuing 
from  blade  A  and  blade  B. 

The  velocities  V^A  and  V*  are  normal  to  the  blade, 
no  account  having  been  taSen  of  the  radial  induced 
velocities. 

The  expressions  for  the  induced  velocities  Vjn  and  Vj 
are  computed  from  a  helix  vortex  having  issued  from  B 
a  point  rH  on  the  blade  (A  or  B)  Z  seconds  before 
(see  Figure  134  ,  formulae  III.1.18  and  III. 1.19). 

The  variables  wZ,  Z  ,  and  are  integration  variables] 
Vf  designates  the  vertical  velocity  of  the  flow  beneatl 
the  disc,  generated  by  the  free  and  connected  vortices, 


Equation  (III. 1.17)  may  then  be  written  t 


( iiij 1 9 20)  n  m 

A 


r  f 

slc^cur0c_li£l:  /  Ha  iS  rt.  <r,r,)olp 

2  di  c  *  dLty  ^  A  '  °  i  2  di  J  Dq  ®  | 


In  theory,  to  resolve  this  equation,  one  may  postulate 

r( r,  „  t  X.  \  <r) 

v  M  ft 

=  'f  <*)  =o 

where  ^(r)  is  a  sequence  of  known  functions  which 
satisfy  the  limit  conditions,  and 

)Tn  are  unknown  constant  coefficients  to  be  determined 

To  this  end,  n  points 


cj  *1, 


are  taken  on  the  blade,  writing  that  equation  111. 1.20 
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is  satisfied  at  those  n  points,  which  in  turn 
gives  n  linear  relations  containing  X,  . .  Y  as 

,  *  i  til 

unknowns.  » 

• 

«  O 

-  This  method  would  call  for  the  calculation  of  2  n*  double 

integrals  (tu£and  ). 

3  P 

These  are  the  integrals  of  (III.1.18),  in  which  -tj: 
is  replaced  by  n  4*  ( )  functions.  These  integrals 
contain  a  fixed  parameter  V^T  such  that  (*fT/R  )%  - 
47X1  /2  PA  Ul  and  n  oarameters  r  =  r«with 

(J  -  1  — n).  '  J 

A  good  approximation  with  this  method  gives  n  equal 
to  5  or  6,  which  in  turn  leads  to  research  for  an 
approximating  method. 

3.12  -  Approximate  expressions  for  the  cores  Na  and  % 

The  core  Ify  reveals  the  discontinuity  for 

r  *  n 


(III. 1.21) 


To  isolate  this  discontinuity,  the  first  half  helixes 
of  the  vortex  O^UiC^ir  must  be  separated  from  the 
core  N^. 


"a  «  ["A.] 


OjZ  s  7T 

<jjZ> 

r 

/ 

'a 

CoZ  s  0 

J  uiZ> 

( II I 4 1.22) 


Let  us  replace  these  first  half-helixes  by  semi-circles 
(  0  and  Vp  *0  )  located  in  the  plane  of  the 

swept  disc. 

The  expression  for  the  velocity  induced  by  a  vortex  cir¬ 
cle  has  been  given  by  Lamb.  Castles  and  Deleew  have  tabu 
lated  this  velocity  numerically. 

In  connection  with  the  method  under  discussion,  the 
following  approximate  numerical  expression  has  been  es¬ 
tablished  for  a  vortex  semi-circle  : 

U  T,  1  _  JLL  _  —h _ 1 

A.  5.7  [  r-<\  r-r*,  r*+*.3Sr* 

Here,  the  first  term  predominates.  The  first  two  terms 
can  be  interpreted  as  two  half-straight-line  vohtricos,  i.e 


/ 
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/"  "ip  A-67 

/  *L  (rt>  rur,t),d.rt  _  n  (rlXl)  rn  (rMi 

Jr.  5  a 


(III 


/ 

.1.25)  J 


If  (r,)  fl  (r.qjdr,  ,  _h  (r,x4)  fM  frMj 

Dr 


*,  *  r.  +  £,  ( rM  -  ro  >  °  <  £-f  <  1 


x2  .  rH  +  t2  (K-  r„  ) 


o  <  <  < 


The  equations  require  elucidation  of  the  following  three 
unknowns : 

1)  -  PH  :  value  of  the  maximum  circulation  intensity 

2)  "  X,  (6,) 

3 )  -  X2  ( £2  ) 

Knowledge  of  rM  is  not  necessary. 

and  rM  can  only  be  estimated,  so  that  one  must 
proceed  by  successive  approximations. 

Though  no  proof  can  be  offered,  it  will  be  assumed  that 
the  abscissae  XA  and  correspond  to  the  marginal 
vortices  observed  in  the  wind  tunnel, 

-  Integrating  with  respect  to  C uZ.  „The  expressions  for 
the  cores  11'^  and  Ng  contain  terms  In  Vp  Z  and  in  cos 
(jjZ  ,  sin  tjj  Z  ,  thereby  excluding  the  use  of  the 
usual  simple  functions.  Vjp  Z  and  c uZ  will  be  treated 
as  two  independent  variably. 

The  first  helix  loop  of  the  vortex,  issuing  from  B 
(0  <  ujZ  4  2W)  (  see  III. 1.18)  will  be  replaced  by  a 

circle  in  a  horizontal  plane  located  beneath  the  blade 
A,  at  a  distance  equal  to  Vf,  Th  .  The  first 
helix  loop  of  the  vortex  issuing  from  A,  i.e. 

7T  4  wZ  4  sir 

will  be  replaced  by  a  circle  in  a  horizontal  plane  loca*< 
ted  beneath  the  blade  A,  at  a  distance  equal  to  VfT  . 
These  vortex  circles  as  p  whole  will  be  termed  "recent 
vortices".  The  velocity  induced  by  these  oircles  is 
given  in  the  tables  compiled  by  Castles  and  De  Leew. 
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The  other  helix  vortices  j 

Vortices  A  3  7T  ^  r 

Vortices  B  2  TT  u;  Z — »< 
;.r 1 1  be-  treated  as  solenoids. 


(III. 1.26) 


&  P 

be  replaced  by  t—  -t —  4z.  (see  equation  III.1.18) 
by  postulating  r*  "‘t 


A  A-and  a*  .  acvft) 

as  vfT  jr,  v 

<v  n 

This  amounts  to  saying  that  a  helix  (or  circle)  JL— 
has  been  distributed  in  uniform  fashion  over  a  0*1 
cylindrical  surface  of  height  Vp.T  (  1/2  \ /e.  T on  either 
side  of  the  initial  helix). 

The  integration  limits  are  then  as  follows  » 

for  U)Z  :  0  and  2  7T 

for  z  -  (giv^n  in  table  below), 

TABLE  3.18 


■i  jr 


Integration  is  now  possible  and  has  already  been  perfor¬ 
med  by  Callaghan  ana  :aslen.  Equations  10,  11  and  12 
in  their  study,  transposed  into  the  annotations  used 

here,  give  the  following  I  ,,  ~ 

vp  C4 

:.1.27)/’_£’M(rr.|drs-/2r  id-  K(t)  (1P,<p] 

J  )rt  "  r’  •  «  J  )'  *VfT  [TT/7ir  k-HVf  8T|  *  J 

vfh. 
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(III. 1.28) 


(III. 1.29) 


(ill. 1.30) 


where  K  is  a  complete  elliptical  integral  of  the 
first  species,  the  modulus  of  which  is  : 

0  2  4  rOi 


(t£  Z)z  +  (r+r4)i 


The  approximate  value  of  K  for  k  values  in  the  region 
of  unity  is  : 

K  (*;  „  Ln  4L  (P.i-**) 


■n  jr 


(Tables  are  available). 

A.  (')>,<)  is  the  Heuman  function  (Byrd's  table). 

For  k  =  0,  this  reduces  to  sin  ,  and  for  k  =  1 

to  2  u> 

—  T 

The  argument  f  is  given  by  : 


tin?  - 


II 

Ci-P 


The  equation  for  the  circulation  T  (III.I.I7)  then 


becomes- : 


/Afp,aTc^Ttvrec"  Jc-fa  f^  +  Vc»+V5A  +^33] 

A 

_  -Lci£t  f  1£  Ld£>  1  f _ 1 _ i  *r<  1 

z  di  J  [fuj[  r.r4  r+r4  ^+2.3 5rf 

ro  J 

where  l/c*  and  are  velocities  induced  by  the  four 
recent  vortex  circles  issuing  from  blades  A  and  B  (oute 
and  inner),  andVi*  and  are  the  velocities  induce 

by  the  four  solenoids. 

Let  rMC  (r4,x,z) 

designate  the  velocity  induced  at  a  point  (  f)  Z  s  0)  on 
the  blade  A  by  a  vortex  circle  of  intensity  PM  . 

The  radius  of  the  vortex  circle  is  JC  (  ac=  X4  and  Xg) ,  lo 
cated  at  a  distance  z  beneath  the  swept  disc  such  that 


*-vf£  f.-r. 


V.T  f, 


r  2.  B  f  r0P  '  A 

The  function  C  is  equal  to  the  quotient  of  a  function 
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tabulated  by  Castles  and  De  Leew  I-75— '  divided  by  the 
radius  X  (  X^and  x2 )  of  the  vortex  circle. 

The  table  below  shows  the  difference  between  the  anno¬ 
tations  used  here  and  those  used  by  Castles  and  De  Leewi 


TABLE  3 .1$ 


Designation 

Present  Text 

Radius  of  vortex  circle 

X  s  X,  or  x2 

Vertical  distance  of  a  point 
(P)  from  vortex  circle 

2-  VtfT 

Horizontal  distance  of  a 
joint  (P)  from  vortex  circle 

r 

Velocity  induced  at  a  point 
(P)  by  a  vortex  circle  of 
intens ity 

r„  c 

Induced  velicity  referred  to 
vortex  intensity 

C(f',3C)X) 

Formula  of  Castles  and 

De  Leew 

xC  C^x.z,) 

Castles  and 
De  Leew 


1.31) 


In  this  case,  the  velocities  and  V/Ci  are  expressed 
by  j 

VCA  s  [C  (r)  XZ  >  Vf  'Z  ^  ^  Xj  »  J 

VCB  *  [C  ^P»  "  C  Cr»  1  Vf  T>] 

Let  Soe  and  Sj  be  taken  to  designate  the  velocities 
induced  at  a  point  on  the  blade  A  by  a  solenoid  of  radiu 
x,  the  density  of  the  vortex  intensity  of  which  is  : 


extends  from  the  swept  disc  to  upstream  infinity 
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(III. 1.32) 


(III. 1.33) 


S j  extends  downwards  from  the  swept  disc  over  a 

distance  z  . 

€> 

In  the  case  of  the  infinite  solenoid,  we  have  : 

1  «  8  *• 

5  Hif  ;<x 

\  0  if  V  >  x 

In  tjje  case  of  the  finite  solenoid,  the  Callaghan  and 
Maslen  result  is  obtained  : 


4  oc  r 


5,(r,X,T  VfT,_< 

Z  'It  4  TTfTSr  4  (x-r;5vfT  1 

h,nfa  — 

’•'hence,  the  velocities  Vg^  and  Vgg  are  given  by  : 


3.13  - 


^  [[4.<hV-  5  (r.xJ.'/fT)] 

-[•S-  <r,*J  ~S  <-r,x4  ■v'Fr,]j 


^SB  [[  X<  *  ■  ^  *  r)x2.  I  t 

f  -|  (r,  x, ) .  3  ( r,  x< ,  Vf  T)J 

The  functions  K  and  \Q  (VJ  k)  exist  in  tabulated  form. 
Estimation  of  Vp  ,  ,  rH  ,  X<  and 

In  the  visualization  tests,  the  basic  length  used  was 

Vf  T 

T  =  2  7T  is  the  period  of  rotation,  and 

~uJ~ 

Vf  the  mean  induced  velocity,  given  by  s 


V 


2pA 


Pj;  was  measured  during  the  tests  and  was  therefore  known 
at  all  times, 

/. 
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In  applying  the  calculation  method,  recourse  may  be 
had  to  the  momentum  theorem: 

(ill. 1.55)  4  7T r  Vp ^  (u>r)  £^(e__YL) 


on  the  assumption  that  Vf  is  being  calculated  for 
r  =  0.5  R  . 

in  mimii  i  --m—  * 

Estimation  of  (rM)t  The  position  of  the  circulation 

m  aximum  must  be  calculated  by'  successive  approximations 
when  applying  the  method.  To  this  end,  a  first  approxi¬ 
mation  is  necessary,  and  it  will-  suffice  to  take  the 
circulation  at  r  =  0.9  R  on  the  basis  of  the  momentum 
equation. 

A  knowledge  of  r*M  is  not  essential. 


Estimation  of  Xj  .  Beneath  the  blade,  the  outer  marginal 
helix  reveals  a  marked  contraction  which  is  a  function 
of  the  load  on  the  disc.  To  simplify  calculations,  one 
may  take  x  =  R  for  the  outer  circles  and  solenoids 
(A  and  B) . 

Estimation  of  Xi  .  The  inner  circles  and  solenoids 
(A  and  B)  can  be  located  in  the  following  manner  :  in  the 
case  of  the  solenoid,  the  radius  of  the  marginal  inner 
helix  increases  markedly  below  the  disc.  In  view  of 
the  variation  in  the  circulation  (see  Section  2),  it 
would  be  logical  to  estimate  the  radius  as  : 


i.e.,  halfway  between  the  inner  extremity  of  the  blade 
and  the  point  where  circulation  is  at  a  maximum.  For 
simplification  purposes,  the  radius  of  the  inner  solenoid 
is  taken  as  half  that  of  the  blade. 

V.'ith  regard  to  the  inner  circles,  their  radius  is  taken 
as  : 


3.14  -  Application  of  the  proposed  method 

It  is  now  required  to  solve  equation  III. 1.50. 

rA  is  the  unknown  function,  while  is  known  approxi¬ 
mately  only.  The  remaining  values  arc  either  known  or 
estimated  ,  as  stated  in  the  previous  paragraph. 


/ 
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Equation  III. 1.50  may  be  written  in  the  following  form  t 

(III. 1.56  rjr,4c^y^(r()_ia_.[T;<r,+  rMTl(r)]_c.  £u;R 

The  terms  (r)  and  Tg  (r)  are  known. 

For  a  given  Pm0  ,  the  function  P  (^)  must  be  found. 
Having  determined  it.  the  P(r)curve  is  plotted  and  the 
P«,  thus  found  compared  with  the  previous  P*i0  »  thereby 
providing  a  better  approximation  Pm* 

This  process  is  pursued  until  a  is  obtained  which 

is  little  different  from  PM 

n  n-n 

It  is  now  proposed  to  show  now  P(r)  is  found  and  bow 
it  is  then  modified. 

-  Resolution  of  the  equation.  Obteinment-of  P(r) 
Changing  the  variable,  we  have  : 

r„  r0  +  ((-tosfj 

(111. 1.57)  )  <ps0  r=r„ 

1  41  -TT  r 

The  unknown  circulation  P^  takes  the  form  ; 

(  p  =  tco  ^-L  ^  21 

(111.1.58)  )  A  2  0  oU  n.i  a 


V  B  0 


r  =  0 


where  c0  is  the  reference  chord,  such  that  c0=  constant 

The  unknowns  arc-  the  dimensionless  coefficients  ya  . 

Performing  the  change  of  variable  in  equation  111.1,36 
integration  can  be  carried  out  easily  by  using  the  wei: 
known  relation  : 


'coj  7 r 

WSV-fftsV  "  s]n  ^ 


(III.  109) 
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6Cl 


Dividing  both  sides  of  the  equation  by  1  C  ■  - 


we  have  : 


2_  r  n't  +  ±C  JL  Z  Yn  n 

n»i  2  dC  11 A  *'"t  ftrv 

=  Ti  (f)  +  rM  t;  (f )  „  e  (Vj 

n  =  6  was  chosen  to  represent  the  function  P  .  One 
nay  therefore  write  that  the  equation  is  satisfied  for 
the  following  values  of  'f. 


30°,  36°,  45°,  60°.  72°,  108°, 
120°,  135°,  144°  and  150° 


The  next  step  is  to  calculate  the  values  of  the  func¬ 
tions  Tj_  (  ^ )  and  T2  (  'f )  for  the  ten  values  of  the 
variable 

This  gives  10  linear  equations  in  #“n  (  n  =  1  ------  6). 

A  number  of  elementary  transformations  lead  to  a  simple 
system  of  equations,  given  in  table  3.20. 

Should  it  be  desired  to  modify  Prt  ,  the  expressions  T^ 

(  ¥')  and  Tg  (  Y )  do  not  change;  only  E  (*f  )  varies, 
but  resolution  of  the  system  III. 1.41  can  be  accom¬ 
plished  quickly. 

The  major  part  of  the  labor  resides  in  calculating 
Ti  ( )  and  T2  (  f  ) .  Tables  would  undoubtedly 
facilitate  this  work. 

The  expressions  for  the  lift  and  the  aerodynamic  moment 
at  the  blade  root  (a  =  0)  are  given  by  the  following 
expressions  s 

’  ^  ^  =  p  (v  r  P 
dr  "  ’ 

n 

-lc0  (u;  R)  (  u;  zL  Jfni«arvf  (  A+2?0-Cos1>) 

^  C  l\ 

n 

*  _f_C  .IS  Xn  SinrtV  l 

l  dr  Z  1  l  J  n.i*  a 
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r  =  r0  +  2^-  M-c«iV) 
dr  ,  sin  f  d  f 


"1 


• 
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TABLE  3.20 

RESOLUTION  OF  EQUATION  (III. 1.3 9) 

The  unknowns  are  ^  ---- 

1 

1 

1 

1 

1 

1 

1 

1 

%  0.58779 

r  (  0 )  -» 

iv  +  k  |  +  2To.  95106 

fl  +  5  k 

1^(36)+e(14 

L  0. 58779]  3 

L  0.58779J  2 

* 

^  0.95106 

1+  k  "f-  £  0.58779 f 

1  +  3  k 

L  0.95106J  3  L 

O.95 106, 

2 

\  0.866 

1  +  k  1  -  y  .0,866  T 

1  +  5  k 

=E(60)+E(120 

L  0.866  J  5  L 

0.866  J 

■  2 

(III. 1.41) 

&  0.866 

1  +  2  1;  1+  K  0.866  T 

I  +  4  k 

=E(30)"E(150 

L  0.5  J  4  L 

0.5  J 

2 

&  0.866 

”1  +  .  2*  1-  X  0.866  r 

1  +  4  k  ~ 

E  +E 

=  (60)  (120 

0.866  J  4  {_ 

0.366  J 

2 

£ 

[i  +  2  k  1-  irf  [i  + 

E 

6  k  1  : 

-E 

_(A5)  (135) 

Z 

L  0.707  J  L 

0.866  J 

2 

k  ,  sp  .  =  i  0  d  C0  . 

r 

IT 

' 7  '  2  di 

R  "  ro 

rnr 

(°)  Uotc  : 

Fr.'-  g  :  ■'  *)  1 ;  *? 5-  •  cl  : r- -t !  • « >,1 3  ,  1 

by  1  +  1  c  dCt  A 

2  dc  VfT 

* 

must  be  replaced 

• 

It 
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3.15  -  Example  of  numerical  application  * 

3.151-The  general  equation  HI.  1.3*0  will  now  be  verified 

by  comparing  tfie  experimental  results  of  Johft  P.  Rabbott 
(see  references)  to  the  theoretic  data  resulting  from 
this  equation,  rotor  characteristics  being  the  same  in 
both  cases. 


The  general  equation,  as  indicated  before,  is  as  followsa 


The  numerical  calculation  of  the  different  terms  -of 
this  equation  shows  that  the  terras  in  j/Wr,  and 

ir<  may  be  neglected  being  indeed  ten  times 
I*1  ♦  r,1  smaller  than  the  term  in  when  oonnec 

ted  whith  sin  \p  and  hardly  reaching  1  of  its  value 
when  connected  with  sin  3  if  ,4^  -  -  10000  -  -  6  . 

Equation  III. 1.30  nay  then  be  written  as  follows  : 


^  (r)  "  je  °  r®c  -  (VCA  +  +  VSA  + 

ft 

aP  1 

3  r,  M-4  '  r-rc 


Its  resolution  will  be  made  with  the  help  of  equations 

III. 1.41. 

J.  F.abbott’s  curve  (T)  was  transformed  in  a 

curve  1  =  f  (  if  )  ,  yj  being  our  variable  in  the  P  equa¬ 

tion. 

The  expressions  permitting  this  transformation  are  : 


and 


with  r  =  0.18, 

O 


JL 

77 


=  f 


u?  r 


F  -•  F0  t  Jzls  (1  -  <&  y) 


A 


The  transformation  of  j.  Rabbott’s  curve  in  curve 
P  =  f  (  f )  »  which  will  be  our  reference  of  comparison, 
is  given  in  Table  3.22  .  See  also  Pig. 135* 

3.152  -  Equation  resolution  hypotheses 

The  characteristics  of  the  rotor  experimented  by 
J,  Rabbott  are  given  in  Table  3.21. 

These  same  characteristics  will  be  introduced  in  the 
considered  equations. 

Plow  velocity  Vf ,  the  knowledge  of  which  is  necessary 
for  computation  of  ,  VqB,  an(*  ^SE*  ls  Given 

by  the  momentum  equation  Vf=  f  \r)  *  at  r  =  0.5 
(fig.  156). 

The  vortices  are  located,  as  indicated  in  paragraph 
3.14,  as  follows  : 

Xz  =  R 

X4  =  0.36  R  for  VCA  and  VCB  and 

0C4  =  0.5  R  for  VgA  and  VgB  . 

For  a  first  approximation  it  is  necessary  to  know  an 
order  of  magnitude  of  PM  , 

In  the  considered  case  PM  =  11  sq.m/sec  will  be  selec¬ 
ted  ,  this  being  the  value  of  the  momentum  equation 
for  r  =  0.9  (see  FiG.136). 

3.153  -  Equation  computation 

Using  Castles  and  De  Leew  Tables  for  computing  V_ 
and  VCB  and  Byrd  Tables  for  computing  VSA  and  VgB 
it  is  possible  to  determine  the  second  member  of  the 
general  equation  III. 1.41. 

Values  of  r  0^  ,  — = —  ,  ■  <  ■  m  r- 

'H  ’ M  'M 

are  given  in  Table  3.23. 

These  values  make  it  possible  to  determine  E  ( Y  ) 
and  E  (  nr  -  f  )  such  as  : 


with 


E  t4  rfj 

®c  and 


-r  -  J_  (  VjA  +V33  »  VCA  +  VC3) 
~  tuK  F\. 


Characteristics  of  rotor  experimented  by  J.Rabbott 
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Item 

Dimension  or  magnitude 

Blade  radius 

R  =  2.30  m  =  7.5  ft. 

Blade  tip  speed 

CUR  =  151  m/sec  =  496  ft/sec. 

Rotor  angular 
velocity 

cu  =  65.5  rad/sec. 

Rotor  Solidity 

6  =  O.097 

Blade  chord 

C  =  0.35  m  =  14  inohes 

Thrust  coefficient 

cT  =  0.00518 

Blade  coefficient 
pitch  angle 

ec  =  9°2  «  0.161  rad. 

TABLE  3.22 


mm 

r 

<14.  lb/.  , 
dr  /loch 

'dr 

P  «9  m  /sac 

S’dagm* 

0.18 

0 

0 

0 

0 

0.30 

0.9 

16.1 

2.84 

45 

0.56 

2.6 

46.5 

4.40 

85 

0.75 

4.8 

86  : 

6.10 

115 

0.85 

7 

125 

7.80 

129 

0.90 

8.1 

145 

8.50 

159 

0.95 

7.1 

127 

7*1 

151 

1.00 

0 

0 

0 

180 

/. 
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Each  term  such  as  vsa  --etc -  is  decomposed 

into  inner  and  outer  VSA,  say,  along  x^  and  xg  . 

In  Table  3.24  are  given  the  values  of  equations  seoond 
members,  that  is  E  (  )  +  E  J [H’-'f  )  and  in  Table  3*25 

the  derived  equations,  in  a  form  still  admitting  P* 
as  parameter. 


3.154 


The  k  coefficient  of  III.1.41  equations  is  equal  to 
0.287 


c  1  d  Cu  TT 

TUr*  “dl"  ha 


0.287 


In  figures  137  and  138  are  indicated  the  variations 
of  VCA»  VCB'  VSA  and  VSB*  versus  . 


In  Table  3.26  are  given  the  values  of  - V©  as 

a  function  of  . 

dc  ns6> 

Determining  r=  lc  tuK  X  ^  sin.  rt'P 


It  is  only  necessary  to  replace,  in  the  equations  of 
Table  3.26,  PM  as  given  by  the  momentum  equation 
for  r  =  0.9  ,  that  is  -:  =  11  sq.m/sec. 

This  first  approximation  gives  ifi - and  lead 

to  the  following  equation  : 


(III. 1.42)  T,  .  =  4.19  sin*f  -  3.81  sin  2  ‘f  +  1.835  sin  3^ 

(sq.m/sec) 

-  0.219  sin  4<K  +  0.118  sin  5  f  -  0.191 

sin  6 

Starting  from  this  equation  we  can  draw  the  curve 
P=  F  (  )  .  Doe  Fig. 139. 

The  Pm  of  this  curve  does  not  correspond  to  the 
selected  PM  and  it  is  thus  necessary  to  begin  again 
the  calculations  with  a  new  Pm  ,  taken  as  a  mean 
value  between  the  two  previous  values,  and  thi*  as  long 
as  the  PM  introduced  into  the  calculations  and  the  P h 
given  by  the  derived  curve  are  not  practically  identica. 

The  convergence  is  rapid  as  shown  by  Table  3.27  where 
are  indicated  the  various  results  of  the  computing 
process. 

Fig. 139  gives  the  variations  of  P=  f  Of )  f°r  the 
different  approximations. 


The  equation  gifing  the  circulation  is  thus  finally 
obtained  : 


TABLE 


U  U 

<d  <y 


s-  $ 

-o 


- T~" 

1 

in 

in 

in 

in 

in 

in 

1 

m 

1 0 

O 

00 

VO 

H 

ON 

00 

CM 

Lf\ 

1 

cm 

CM 

m 

r- 

s 

00 

ov 

o\ 

O  1 

1 

t 

1 

• 

0 

• 

O 

• 

0 

• 

• 

0 

• 

0 

« 

0 

• 

0 

• 

0 

• 

0 

150  0.153  0.216  0.0166  0.119  0.0094  0.260  -  0.033  0.21  -  0.0396 
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TABLE  3.24 

% 

E  <<f>) 

30 

• 

E  (30)  =  0.0378  -  0.000837  PM 

56 

E  (36)  =  0.0419  -  0.00181  Th 

45 

E  (45)  =  0.0484  -  0.00252  rM 

60 

E  (60)  =  0.062  -  0.00474  rM 

72 

e  (72)  «  0.075  -  0.00627  r„ 

108 

E(103)  =  0.115  -  0.00725 

120 

E ( 120)  =  0.128  -  0.00702  PM 

135 

E ( 13 5 )  =  0.142  -  0.00646  rM 

144 

E ( 144 )  =  0.149  -  0.0057  rM 

150 

e ( 150)  -  0.153  -  0.00501  rM 

TABLE  3.25 

System  of  equations  £o  be  solved. 


with  r  as  parameter 
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(1)  1.86  +  2.85  *u 

(2)  1.441 ^  -  2.01 

(3)  1.812  y2 

(1) '  1.151*, 

(2) '  0.874^  +  2.}45  y3 

(3) '  1.24  fc,  -  1.12  *3 


=  -  0.0576  +  0.002086 
=  -  0.033  +  0.00114  Prt 
-  3.44*6  =  "  0*0468  +  0.00197  Pm 

.  2.3*  =  0.095  -  0.00588  Pm 

=  0.09502  -  0.003755  Prt 

=  0.095  -  0.00676  Ph 


TABLE  3.26 

Determining  ^  versus  P m 


*n  are  dimensionless,  PM  being 

Y,=  0.0849 

~ 

0.00515  Pm 

v2=  -0.0268 

+ 

0.00095  Tm 

0.00891 

+ 

0.000318  Trt 

*4=  -0.0027 

+ 

0.000111  Prt 

*5=  0.00113 

- 

0.0000304  rM 

*e=  -0.000495 

- 

0.0000728  PM 

sq.m/sec . ) 


*.84 


’  ,  =  6.275  sin  -  2.804  sin  2<P 

(sq.m/sec) 

(III. 1.45)  +  1.708  sin  5  V  -  0.2637  sin  4  +  0.130  sin  5  V  -  0.162 

c 

sin  6  f 

It  is  now  possible  to  compare  this  curve  r.((V)  , 
as  given  by  equation  III. 1.43.  with  the  curve  obtained 
by  J.Rabbott.  (See  Pig. 135  and  140). 

In  the  same  way  as  previously  the  curve  Z"1  =  f(if)  may 
be  transformed  backwards  into’  '=  f  (?)  and  compared 
to  the  original  J.  Rabbott’s  curve.  See  Fig.l4l. 

It  may  be  seen  that  the  comparison  of  the  two  curves  is 
quite  satisfactory. 

In  Pig. 142  is  given  the  final  curve  P=  P( ¥*)  . 

This  Fig.  shows  also  all  the  components  of  P  as  indica¬ 
ted  by  equation  III. 1.30. 

3.16  -  Attempts  to  use  simplified  methods 

In  this  section  are  analysed  simplified  methods,  the 
objective  of  which  is  to  reduce  the  number  of  terms 
of  the  equation  giving  the  circulation. 

First  simplified  method. 


I 


(III. 1.44) 


The  first  two  terms  of  the  second  member  of  this  equa¬ 
tion  are  the  same  as  previously. 

The  last  integral  means  that  to  the  helix  vortices  loca¬ 
ted  under  the  swept  disc  has  been  substituted  an  infinity 
of  elemental  solenoids  (A  and  B)  whoso  number  of  turns 
per  unit  coil  length  is  : 

_i_ 

vfT  ac, 

It  may  be  noted  that  the  lower  limit  of  the  last  Integra] 

is  r  instead  of  r  because  inside  the  solenoid  the 

0  . 


001185  -  0.175 


001097  -  0.162 
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axial  magnetic  field  is  constant  and  it  is  null  outside 
the  solenoid. 

Let  us  perform  the  transformation  already  indicated  t 


r  =  rc  +  *~r°  (^-cos'f ) 

z 

n 

r  a  21  A*  $irv 


The  above  equation  may  then  be  written  as  follows  s 


(III 


h 

H  A  Sin.  nV  ^f4.c  4$  _±  +  i  ci$r 
.1.45)  r*s‘1  *  dt  VpT  *  d.1  R_£  -H.4  sia^P' 

-  Tc  $r<*\rec 


The  system  of  equations  III. 1.24  is  still  valid  if  1 
is  replaced  by  (  1  +  1  c  .  d  d  4  *  and  if 

./*  2  "ST  '  vpT 

E  ('f>  “  ic 

For  the  numerical  application  selected  previously 
(j.Rabbot's  rotor)  and  using  the  metric  system,  the  folio 
wing  results  are  obtained  : 


ic  =  2- 

A  c  A  Ql  L  =  1.04 

Z  di.  R_  r0 

Ic  -  2  .  IE  «*  0.287 

2  <1^  R.re  11.4 


=  23.75  (  ec  =  9*2) 


Vf.T=  0 .74- 


1  +  I,  -±_ 

2  °  die  Ve  T 


»  2.325 


The  values  of  E  (  ¥  )  =  1  c  A  *•<-  cor  0C  are  given  by  Tabl 
3 1 28  2  di. 
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TABLE  J.  28 


£  (>P)  =  £c  un0 1 


36 

45 

60 

72 

108 

120 

135 

6.125 

7.12 

9.15 

14.5 

16.85 

■ 

18.9 

20.9 

The  resolution  of  the  system  of  equations  gives  : 
r  ( 'P )  =  6.38  s  in  ^  -  2.34  sin  2^+  O.99  sin  3  f 

-  0.38  sin  4*f  +  0.189  sin  5  V  -  0.0812  sin  6  ¥ 

The  harmonic  analysis  'Runge's  method)  of  J.Rabbott’s 
experimental  curve  gives  : 

P(^)  =  6.05  sin  4*  -  2.46  sin  2  +  1.56  sin  3  f 

-  0.511  sin  4  H5  +  0.192  sin  5  ¥  . 

Fairly  large  discrepancies  may  be  observed  for  the  terras 
in  sin  4*  ,  sin  3  'f*  and  sin  4  though  the  sum  of  terms 
is  practically  the  same  :  7.6I  against  7.59. 

The. corresponding  curve  is  given  in  Fig. 143. 

Second  simplified  me t h od 


PA=  1  c  i£t  wr  et  _  ±  c  _i£fc  _L 

t  d-C  !  Ili  41T 


(III. 1.46) 


_  ±  C  Vr 

Z  di  r 


4  TTC  vf 


co  r 


k  c  Ag>-  Ctur  0,  _  Vr 


The  second  expression  is  the  momentum  equation. 


The  successive  transformations  show  a  relationship  with  th 
equation  of  che  infinite  solenoid. 


r  _  r. 

A  B 


DC 


cU 


4  7Tr  W 


p  Vp  =  co  r  (  ^  +  Tj  ) 


cu  =  UL 

T 


Vp  = 


Ztt  r 


4  tt  r  Vp  T 


+  i) 


Vc  = 


2vfr 


(  PA  +P 


B  VCT  A 


V/> 


or  _j_  d. 


3r  vfT 


The  change  of  variable,  from  r  to  f  ,  leads  to  the 
following  equation  s 


re 


T  s  X  Aa  sin  a'P 

rt»* 


a 


(111.1.47) 


n«4 


V"v 


in  ,  ,  4  d  Cl  ^  rv  . 

M  +  ITT  35T  ; 


=  'ci^wre,  _4 

2  dC  c  2 


±c  V 

*  '  4l 


The  III. 1,41  equations  system  remains  valid  if  a  substitu¬ 
tion  is  made  for  s 

£  a  -1  c  cur  ec  _  4  C  iii  7C 

2  die  c  T  <U  F 


Taking  again  the  previous  numerical  application  ,  the  foil 
wing  is  obtained  (metric  system)  » 

/. 


1  c 

2 


1 

T 


0.245  # 
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1 

2 

1 

2 


Ml 

di 


C,  uj  “R  r  8^ 

dLo 

c  ^0  =  0.98 

dLC 


=  25.75  r  (  ec  =  9?2) 

O 


TABLE  5 . 29 


a 

56 

45 

60 

72 

108 

120 

155 

144 

150 

5.59 

6.14 

7.14 

9.15 

14.5 

16.9 

18.9 

m 

mm 

■mi 

mm 

4.61 

4.95 

5.5 

6.5 

8.7 

9.6 

10.5 

11.1 

11.5 

_ 

11.5 

_ 

4.52 

4.85 

5.59 

6.57 

8.52 

9.41 

10.2 

10.88 

m 

Bl 

I.07 

1.29 

5.98 

7.49 

8.70 

10.02 

10.95 

11.69 

Resolution  of  type  III. 1,41  system  of  equations  gives 
the  following  results  : 

r  (H>)  =  6.05  sin  i  -  2.62  sin  2  +  O.505  sin  5  f 

-  0.52  sin  4  +  0.469  sin  4  V  +  0.501  sin  5 


The  corresponding  curve  is  given  in  Fig. 145. 

In  this  Figure  are  grouped  for  comparison  s 

-  the  experimental  J.Rabbott’s  curve, 

-  the  curve  derived  from  the  complete  theoretical 
method, 

-  the  two  curves  derived  from  the  simplified  methods. 

The  proposed  complete  theoretical  method  shows  to  be  far 
the  best.  It  may  be  used  to  determine  the  loads  acting 
on  the  blades. 


The  first  simplified  method  is  less  good  'but  it  may 

be  used,  however,  for  overall  performance,  calculations  * . 
The  proposed  complete  method  maj’  s.ill  De  improved  by 
taking  into  account  the  contraction  of  the  outside  helix. 
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3.2  -  Problem  of  the  non-stationary  conditions 
applied  to  helicopter  rotors. 

© 

Attempts  to  find  an  approximate  solution. 

*  3.21  -  Objective 

The  local  lift  of  a  blade  In  non-statlonary  conditions 
may  be  considered  a  the  sum  of  two  effects  » 

-  circulation  due  to  a  distribution  of  velocities,  normal 
to  the  (thin)  airfoil  section  , 

-  impulsion  of  air  on  the  airfoil  due  to  accelerations 
(derivatives  of  those  normal  velocities). 


Defining  the  circulation  is  difficult  because  of  the 
presence  of  the  free  vortices  escaping  from  the  trailing 
edge  (wake),  the  intensity  of  which  must  satisfy  the 
vortex  conservation  law. 


These  free  vortices  (termed  radial  vortices,  by  opposi¬ 
tion  to  helix  vortices),  after  appearing  at  the  trailing 
edge  of  blade  A,  move  away  from  the  latter  during  a  half 
period  (except  in  the  inversion  circle)  then  come  back  unc 
the  rotor  disc  and  so  forth. 


er 


Those  of  the  radial  vortices,  which  are  formed  less  than 
half  a  period  ago,  will  be  termed  recent  vortices. 

The  totality  of  the  radial  vortioes  escaping  from  the 
blade  trailing  edge  will  then  be  separately  considered 
according  ,/hether  they  are  recent  vortices  or  not. 

The  recent  vortices  will  be  considered  as  the  vortices  of 
the  non-stationary  wake. 

The  other  radial  vortices,  whether  they  come  from  blade  A 
or  from  blade  B  ( in  the  case  of  a  two  bladed  rotor)  will 
be  considered  as  inputs  in  the  same  way  as  the  blade  sec¬ 
tion  pitch  angle  variations. 

In  this  way  it  becomes  possible  to  make  use  of  number  of 
results  of  the  fixed  wing  non-stationary  theory.. 

However,  these  results  lead  to  expressions  which  are  diff :  - 
cult  to  handle  when  applied  to  helicopter  rotors,  and  for 
this  reason  approximate  formulae  will  have  to  be  derived. 


In  the  first  part  of  this  section  will  be  reminded  some 
basics  of  the  non-stationary  theory. 


In  the  second  part  will  be  exposed  an  approximate  method 
used  for  fixed  wings. 


/ 
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3.22 


In  the  last  part  these  results  will  be  applied  to  helicopt 
rotor  blades. 


elr 


Basics 


3.221  -  It  is  not  our  intention  to  expose  here  the  non-stationa- 
ry  theory  of  lifting  wings'but  merely  to  give  some 
results  which  might  be  useful  for  a  better  comprehension 
of  the  following  text. 


(III. 2. 48) 


The  lift  of  a  blade  section  dr,  either  in  stationary  or 
in  non-stationary  conditions,  is  due  to  a  pressure  dif¬ 
ference  between  the  upper  and  lower  surfaces  of  the  air 
foil  section.  It  may  thus  be  obtained  by  integrating 
the  pressure  differences  along  the  airfoil  chord.  We 
may  write  : 


dL 

dr 


)  dx 


p^and  designate,  respectively,  the  local  pressures 
on  the  airfoil  lower  and  upper  surfaces  (See  Pig. 144). 


In  Pig. 145  are  given  characteristic  pressure  distri¬ 
butions,  for  stationary  (a)  and  non-stationary  (b) 
working  conditions,  in  the  case  when  the  airfoil  may 
be  considered  as  a  plane  thin  sheet. 


Presence  of  large  "bumps", in  the  rear  part  of  the  air¬ 
foil  section  (rigid  and  flapless)  indicates  rates  of 
non-permanence  of  the  motion. 


Their  absence  indicates  that  the  working  condition  is 
quasi  stationary. 

3.222  -  The  study  of  non-stationary  conditions  is  made  easier 

by  substituting  to  the  notion  of  section  angle  of  attack 
a  notion  of  distribution  of  the  velocity  Vz  (x,t)  normal 
to  the  airfoil.  Vz  varies  along  the  chord  (x)  and  with 
time  (t). 

The  existence  of  two  typical  distributions  of  Vz  (trans¬ 
lation  and  rotation)  gives  rise  to  a  circulation. 

The  existence  of  accelerations  yy  Vz  (x,t)  gives 
rise  to  an  impulsive  lift,  i.e.  to  the  effect  on  the 
blade  of  the  accelerated  air  s irrounding  it. 

In  stationary  conditions,  lift  is  due  to  the  circulation 
r  exclusively. 


/. 


(III. 2. 49) 


yr= 
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Circulation  is  kept  up  by  the  VzQ  component  of  the 
velocity,  normal  to  the  airfoil  (see  Pig. 146).  This 
velocity  is  constant  in  time  /  5  l/r  ,o)and  along  the 
chord  (^Lv^.o).  1 

All  happens  as  if  the  airfoil  section,  placed  in  an  air 
stream  VR,  was  translated  with  a  speed  Vs0,  relative  to 
the  surrounding  air,  in  a  direction  normal  to  the  ohord, 


(III. 2. 50) 


r  =  _L  C.  V- 

2  cLL  2 


j  r  dCL  V/to  y 

T  TTTf  R 


>/e  will  say  that. the  airfoil  section  placed  in  an  air 
stream  is  a  black  box  whose  input  is  a  translation  velo¬ 
city  Vz0  and  whose  output  is  the  circulation. 

For  a  helicopter  in  hovering  conditions  : 

V«£o  a  Co  C  9  -Vt 

The  velocity  wr  9  is  constant  in  time  and  along  the 
chord. 

The  velocity  V*  induced  by  the  helix  vortices  is  oonstan 
in  time  but  not  absolutely  constant  along  the  ohord.  In 
fact  it  is  the  velocity  of  the  neutral  rear  point  (loca¬ 
ted  at  three  quarter  chord  from  the  leading  edge). 

As  it  has  been  said  before,  in  non-stationary  conditions 
lift  is  produced  by  a  double  effect  : 

-  the  circulation  due  to  the  Vz  (x,t)  velocities, 

-  the  impulsive  force  due  to  accelerations  —  \L. 

dt  z 

Let  us  first  examine  the  circulation  . 

It  is  convenient  to  write  j 


(III. 2. 51) 


Xc  ^  cos  9 


x’-f 

X.*i  4.  iL 

Z 


9  r  7T 


leading  edge 
trailing  edge 


and  to  express  Vz  (x,t)  =  Vz  (  0,t)  in  the  form  of  a 
Fourier  series  (oosine  series). 


(III. 2. 52)  Vj  (0, t)  s  V-  (t)  4*  X  VZfv  (k)cos  r\0 
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The  first  term  of  the  second  member  of  this  equation 
corresponds  to  a  translation  of  the  airfoil  section 
along  the  normal  to  this  section. 


The  second  term  : 


cfc)  cos  e  s  _ 


*4* (b) 


represents  a  rotation  of  the  airfoil  section  around  a 
point  located  at  half  chord. 

All  the  components  of  the  V2  velocity  give  rise  to  a 
pressure  distribution  (Thin  airfoil  theory,  by  Munk- 
Glauert),  but  only  the  two  first  components,  i.e  * 

V^o  <fc;  translation 

(t)  rotation 

give  rise  to  a  circulation  around  the  airfoil  which  gene¬ 
rates  lift. 

For  the  other  terms  of  V2  ,  such  as  n  >  1,  the 
(III. 2. 43)  integral  is  equal  to  zero.  ' 

We  may  note  that  the  translation  velocity  Vz  (t)  is  not 
a  mean  velocity  along  x  : 


VZ0  (k)  ¥=  “f  /  ^  (X, 


k  )  dx 


It  is  the  constant  term  in  the  cosine  series  i 


(t)  -  4 

7T 


'  C/2  J  ^ 

vz  ( x,k >  — 


We  will  admit  that  the  discrepancy  is  small „ between  these 
two  values. 

The  effect  due  to  accelerations  will  be  examined 
later. 

Schematic  description  of  the  circulation.  Relative 
importance  of  translation  and  rotation  motions. 

Let  us  consider  an  airfoil  section  placed  in  a  bi-dlmen- 
sional  flow, with  a  uniform  wind  VR  (t)  which  may  vary 
in  time. 

The  possible  inputs  are  : 

a)  a  translation  of  the  airfoil  section  along  a  direction 

/. 
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normal  to  this  section*:  7  (t). 

O  • 

b)  a  rotation  (t)  around  a  point  located  at  half 
chord. 

The  Vz  (x,t)  velocity  is  then  : 


(III. 2. 53)  (x,t>  =  [j  ft)  +  ji-  [V,  (t)  w]  J  +  «  <t)  a: 

If  no  vortices  existed  at  the  trailing  edge,  the 
corresponding  circulation  would  be  : 


(in.2.54)  r\t>.  -f  C  [j  (t,  +  -3-  (VR«o] 

As  circulation  varies  in  time,  vortices  having  the  same 
intensity  and  their  signs  opposed  to  the  successive  in¬ 
creases  of  circulation  around  the  airfoil  will,  in  fact, 
be  escaping  from  the  trailing  edge,  in  virtue  of  the 
vortex  conservation  law. 

These  new  free  vortices,  parallel  to  the  trailing  edge 
and  orthogonal  to  the  stream-lines,  will  induce  a  new 
distribution  ^Vz(x,t)  around  the  airfoil  which,  in  turn, 
will  modify  the  circulation  and  give  rise  to  a  new  depar¬ 
ture  of  free  vortices,  and  so  forth...... 

If  we  consider  again  our  image  of  the  "black  box",  des¬ 
cribed  previously, the  input  is  the  Vz  given  by  equation 
(III. 2. 55)  buttt5utput,  i.e.  the  circulation,  is  not  P* 
given  by  equation  (III. 2. 54),  but  ,  given  by  j 

(III. 2.55)  H  =  □  P  * 


where  □  is  an  operator  of  P  accounting  for  the  passage 
through  the  black  box.  1 


A  vortice  -  dP  (see  Fig. 148)  which  has  escaped  from 
the  trailing  edge  Z  seconds  before  the  present  time  t, 
is  located  at  a  distance  : 


(j-  -  -r  +  /\< 

Jo 

It  induces  at  a  point  x  of  the  airfoil  a  velocity 

as  : 


sud 


/, 


^  (x,b 


£ 

c 


dr  _ ± _ 

2TT  if-X 


4  -t-CoS© 


By  braking  down  this  expression  into  a  Fourier  harmonic 
series  ,  we  obtain: 


(111.2.56)  VZ(X,\:)  = 


dr 

27T| 


n. 


M) 


i^L^coa  ne 

/I 


and  using  the  transformation  indicated  by  Rufus  Isaacs: 


cosne  de  ,  7r 


a-cos6 


/a- 1 


the  local  vortex  distribution  may  be  written  : 


2b  <e't> 


_d£_ 

27,f 


L 


coUn.JL_4  (-1  jin 


Among  other  terras,  we  can  recognise  in  (III. 2. 56)  the 
presence  of  the  translation  and  rotation  terms. 


As  the  total  circulation  around  the  airfoil  is  only  due 
to  the  translation  and  rotation  terms,  the  modification 
of  the  circulation  around  the  airfoil  by  a  wake  free  vort  t;, 

located  at  a  distance  S»  -£L  f rom  the  half  chord,  is  : 

s  2 


(III. 2.57) 


or  : 


Z7 r  ( 


( 


additional 

This7£irculation,  highly  intensive  at  the  trailing  edge 
(  £  =  1),  becomes  very  rapidly  equal  to  zero,  as  the  vorte 
moves  away  (  £  — *•  00  ) . 


We  may  note  that  it  is  not  only  a  function  of  the  intensi 
ty  of  variation  of  the  circulation  d  P  but  also  a  /. 
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function  of  the  velocity  VR  ,  with  which  the  free  vor¬ 
tex  is  carried  away  by  the  air  stream  (see  Table  5.50). 

No  great  error  would  be  made  if  the  effect  of  the  wake 
was  modified  or  even  suppressed  beyond  a  distance  of  4 
to  5  chord  lengths. 

The  circulation  produced  by  the  translation  effect  is 
always  larger  than  the  one  produced  by  the  rotation. 

The  ratio  between  those  two  effects,  is  as  follows  t 

translation  £[  >A 

6  rotation 

The  values  of  this  ratio,  versus  £  ,  are  given  in 

Table  3.50. 

f;  times  the  half  chord  length  being  the  distance  of 
the  free  vortex  to  the  middle  of  the  chord. 


TABLE  5.30 


4 

1 

2 

4 

5 

10 

/£  - 

0.73 

0.29 

0.22 

0.11 

1 

0.27 

0.13 

0.10 

0.05 

Figure  145  (b)  gives  (ref.  Von  Karman)  the  pressure 
distributions  generated  by  the  wake  vortex,  as  a  function 
of  its  distance  to  the  airfoil. 

We  may  note  here  that  ty*  (t)  intervenes  in  two  ways  s 
on  the  one  hand,  it  appears  in  the  expression  giving 
(see  equation  III. 2. 53),  whence  it  plays  the  role  of  an 
input,  and, on  the  other  hand,  it  appears  in  4  , 

which  means  that  in  the  circulation  black  box  it  plays 
the  role  of  an  impedance  having  an  influence  on  all  the 
Inputs. 

A* 
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3.224  -  Approximate  configuration  proposed  for 
helicopter  rotor  blades. 


In  the  case  of  a  helicopter  rotor  blade  the  free  radial 
vortices  (wake)  draw  a  spiral  shaped  surface.  They 
first  move  away  from  the  blade,  during  half  a  period 
(except  in  the  inversion  circle)  then  come  back  under 
-the  rotor  disc,  in  the  vicinity  of  the  blade. 

In  order  to  make  maximum  use  of  the  results  available 
for  fixed  wings,  the  wake  of  the  blade  will  be  cut  down 
into  two  parts. 

To  the  first  portion  of  the  spiral  shaped  sheet,  which 
was  formed  between  half  a  period  ago  and  the  present 
time  t,  will  be  substituted  a  plane  sheet  moving  away  to 
infinity,  because,  as  it  has  been  said  before,  its  in¬ 
fluence  after  four  or  five  chord  lengths  is  small. 

The  other  portion  of  the  spiral  shaped  sheet,  formed  at 
an  earlier  time,  will  be  regarded  as  an  effect  which 
would  be  foreign  to  the  wake,  comparable  to  a  blade  pitcl 
angle  variation. 

3.225  -  Accelerations  effects 

In  the  case  of  a  two-dimensional  flow  (infinite  aspect 
ratio)  the  lift  generated  by  impulsion  is  equal  to  the 
inertia  force  of  a  mass  of  air  whose  volume  is  a  cylin¬ 
der  having  c  for  diameter  and  the  wing  span  for  heighlj 

The  acceleration  is  : 


(JII.2.58) 


dfc 


3 

i(t'  f  TT 

Z  (t)  +• 

.  0  3oi 

3v„ 

'  W  + 

Jji 

3t* 


In  the  case  of  a  wing  having  a  finite  aspect  ratio  the 
above  force  is  to  be  multiplied  by  a  coefficient  which 
is  the  ratio  between  the  span  and  the  half  perimeter  of 
the  wing. 

In  the  case  of  a  helicopter  rotor  this  ratio  will  be 
taken  equal  to  one  . 

3 . 23  -  Approximate  method  used  for  fixed  wings 
3.231  -  Objective.  Limitative  hypotheses. 

The  wing,  placed  in  a  wind  Vr  (  t)  is  considered  as  an 

- * - — _ A 


* 
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impedance  system  Z  (p)  whose  input  is  the  velocity 
Vj  (x,t)  normal  to  the  airfoil  (except  in  the  zone  of 
the  trailing  edge  wake)  and  whose  output  is  the  circula¬ 
tion. 

This  system  is  generally  governed  by  an  integral  equation 

It  may  be  theoretically  replaced  by  a  systenrof  ah 
infinite  number  of  differential  equations. 

A.  first  approximation  is  thus  made  when  this  system 
is  replaced  by  a  single  differential  equation. 

The  second  approximation  comes  from  the  non  linearity. 

If  the  airfoil  is  stalled  the  theory  is  not  valid  any 
more  and  another  theory  is  necessary. 

As  mentioned  previously,  the  velocity  V/f*  (t)  plays  a 
double  role  in  the  circulation  :  firstly  as  an  input, 
secondly  as  a  wake.  The  influence  of  VR  on  the  ‘mpedance 
is  thus  a  cause  of  non  linearity. 

It  may  be  noted,  however,  that  if  the  local  advance  ratio 
p  is  small  the  error  made  in  considering  the  system 
as  linear  is  also  small. 

We  will,  therefore,  hence  consider  the  system  as  being 
linear. 

In  these  conditions  we  are  led  to  the  following  problem  i 
we  have  to  define  the  impedance  of  a  system,  by  means 
of  a  linear  differential  equation,  knowing  the  latter's 
response  to  : 

-  a  sinusoidal  input 

-  a  unit  step  input. 

3 . 23 2  -  Determining  an  approximate  solution 


Let  us  consider  a  wing  of  aspect  ratio  A,  placed  in  a 
wind  speed  Vp  (t). 

The  input  ^  is  an  unit  step,  i.e.,  being  zero 

at.  time  0  it  suddenly  takes  the  value  V Zo 

then  remains  constant.  V* cfc) 

This  case  has  been  studied  by  Wagner  for  infinite  aspect 
ratios  and  by  Jones  for  aspect  ratios  of  6  to  5. 


The  approximate  unit  solutions  given  by  Jones  are  : 
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['i-'R(xJ)] 

—  0.0 4Si  _  0.3/4 

=  1-0.165  c  -  0. 335  a. 

At  00 

r  ~  i 

-<$4 

[  1-  R  (^1)J 

=  -1  -  R0  ^ 

A 

dLj  s  *¥*A ll 

c 


The  values  of  RQ  and  S  ,  for  A  *  5  and  6,  are 

given  in  Table  3.31 


A 

5 

3 

0.283 

0.540 

6 

0.361 

0.381 

The  approximation  made  can  only  be  estimated  forAroo  ; 
the  Laplace  transform  of  the  response  to  the  unit  step 
input  is  compared  to  Theodorsen's  curve  (see  Fig. 149). 

lie  will  admit  that  this  approximation  is  valid  for  values 
of  A.  such  as  3  and  6. 

For  a  helicopter  rotor  we  will  admit  j\.  =  6.  In  this 
way  we  will  not  only  take  into  account  the  vortex  paral¬ 
lel  to  the  trailing  edge  but  also  portions  of  lateral 
helixes. 

Formula  III. 2. 5  for  A  =  6,  may  be  construed  as  the 
solution  of  a  differential  equation.  We  obtain  i 


(III. 2. 60) 


clt  * 

&A 


( -f  -  )  -+*  <5  i  * 


1  *  is  the  static  angle  of  attack  (similar  to  I"1*) 
i  is  the  dynamic  angle  of  attack  (similar  to  P  ) 


/. 
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It  is  to  be  noted  that  we  may  substitute  to  1,  either 
C,  =  or  JP__ 

L  <xi  VR 

The  equivalent  expression  may  then  be  written  as  follows i 
£  0^  .  (J  C  dl  d.  I  p  .  S  •  X 

is  dT  +  =  IT  +Sl 

—  —  )  +  <5r-c  i  ,  P* ,  , j  „  >  jn 

«  it  (S,,+  -  T  w  l~vT }  1  *')  +dr 

K 

f  *  '*4*  -  £)"•«-<.,  f.  (4s-(.-i5,)4jr' 

=  -1- VR  tr,t) 


If  the  derivative  is  considered  as  negligible  and 

Vr  considered  as  constant  during  a  range  of. four  chord 
lengths,  it  may  be  said  that  the  wake  vortices  move  with 
a  constant  speed. 

This  result  is  only  valid  for  the  translation  motions. 

In  order  to  introduce  the  rotation  we  will  calculate 
an  equivalent  translation  velocity  4  Vr  by  means  of 
the  following  formula  i  0 


(III. 2. 62) 


AV  s  u  (_£._  £L  ) 
z  V  4  &  ' 


a.c  defines  the  position  of  the  center  of  rotation, 
taken  as  positive  in  the  direction  starting  from  the 
middle  chord  point  towards  the  trailing  edge.  ‘ 

If  rotation  takes  place  around  the  middle  chord  point 
then  a  =  0. 

If  rotation  takes  place  around  the  forward  neutral  point 
(aerodynamic  center)  we  obtain  t  a  =  -  l  and  a  o  »  -  c 

2  ~  X  * 

5 • 24  -  Application  to  a  two-blade  helioopter  rotor 
5.241  -  Disturbed  working  conditions 
5.2411  -  Definition 

Let  us  consider  that  hovering  conditions  is  the  initial 
state  of  operation  of  the  rotor  and  that  small  dis- 

_  /. 
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(III. 2. 63) 


turbances  are  then  introduced. 

These  disturbances  may  be  classified  in  two  types: 

-  symmetric 
*  -  antr‘ symmetric 

(see  Tables  3*32  and  3.33). 

The  symmetric  Derturbances  increase  rotor  lift  and 
modify  the  rotor  conicity  angle  a0. 

The  antisymmetric  perturbances  do  not  modify  rotor 
lift  but  tilt  the  rotor 

’  fbit)  =  <30  (fc)-^(t)  co$  Y- (t)3in  Y . 

I.  J.  M.+  M  1  cos  Y  +  sin 

Jlis  +  cu  Jh  i!£f_  +0O*A0  =  (  t) 

dt*  I.  db  °  I  ; 

i'  +  CU  .-Si  ^  tiw  b,  JLi-  c  --^1  (t) 

r.  x#  i, 

2c +■  co^  £i-  _  CD  ^3_  =  (t) 

^  5  *r  f  ijt  rridr 

TABLE  3.32 

Symmetric  Disturbances 


Disturbance  originated 
by  : 

Variable 

Disturbed 

state 

Collective  pitch 

ec 

6C  +  A6c 

Rotational  speed 

to 

CU  4-  A  CU 

Ascensional  velocity(#) 

Vi 

e  +*vt 

(°)  not  experimented. 
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TABLE  5.: 


Antisymmetric  -  Disturbances 


Disturbance 
originated  bs 


4  Cyclic  pitch 


5  Horizontal  speed 


6  Rotor  tilt  or  angle 
of  attack 


Variable  Disturbed 
state 


sin  V  0  +  A$z 


oi  (oL) 


Every  disturbance  modifies  the  lift  of  a  blade. 

Let  us  assume  that  this  modification  consists  in  an 
increase  of  lift  of  blade  A. 

The  increase  which  appears  at  the  birth  of  the  dis¬ 
turbance  is  attenuated  by  the  blade’s  own  recent  wake 
effect  and  by  the  impulsion  (Z0i) .  As  lift  increases 
the  hinged  blade  comes  up  (£  )  and  this  motion  gives 

rise  to  a  double  effect  : 

a)  the  flapping  velocity  opposes  lift, 

b)  as  the  blade  comes  up  it  gets  away  from  the 
former  helix  vortices. 

This  decreases  their  influence  on  the  blade 
and  ,  therefore,  decreases  blade's  lift. 

After  approximately  half  a  period  the  Influence 

of  the  radial  vortices  issued  from  blade  B  becomes 
effective  on  blade  A,  first  giving  rise  to  a  deorease 
then  to  an  increase  of  lift. 

Finally  intervene  the  marginal  vortices. 

The  latter  do  not  take  their  position  immediately  but 
the  different  helix  portions  get  into  place  progressi¬ 
vely  as  they  are  being  reached  by  the  disturbances  or 
the  vortloes. 


A  movie  film  has  been  taken  of  this  phenomenon  (see 
reference  8)»whioh  illustrates  what  has  been  said  above 

Let  us  compare  now  the  time  constants  Zo  , 

Account  being  taken  of  the  wake  and  of  the  impulsion,  / 
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the  lift  of  blade  A  becomes  established  when  the 
blade  has  travelled  over  a  distanoe  of  n  chord  lengths 
(say  ,  n  =  4) . 


(III. 


(III. 2 


2.64) 


TCi  =  .  UL  n.  bc_  .  J_  jl  (?.-£-) 

'  tur  ai  TTR  P  zb  1  TT1 


For  a  two  tladed  rotor,  having  a  solidity  Of  =  0.10, 
the  time  constant  is  : 


that  is  s 


(  T.  J2) 

U) 


0.10  T  for  the  blade  tip, 

0.153  T  for  the  section  at  O.75  R 

0.20  T  for  the  section  at  0.50  R 

and  T  for  the  section  at  0.10  R. 

In  the  case  when  the  influence  of  lift  on  the  flapping 
motion  is  to  be  calculated  £t  may  be  taken  that  the 
non  stationary  effect  is  the^which  concerns  the  section 
at  0.75  R. 

V/ hen  it  is  wanted  to  study  lift  distribution  along 
the  blade,  the  variation  of  the  non  stationary  effect 
along  the  blade  radius  may  be  taken  into  account. 

The  flapping  time  constant  Z  ( 02)  is  as  follows  (see 
equation  III. 2. 63): 


,65) 


JJL  rJL  )c  T  W  I. 

w  1  *  j  Tr  gJ- 


For  helicopter  rotor  blades  the  ordinary  value  of 
is  comprised  between  0.5  and  1.25.  I* 


Whence  j 

0.8  4  Gj  4  2.0 

The  value  of  Z0 ^  is  therefore  comprised  between  a 
quarter  and  two  third  of  a  period  (we  leave  cut  the 
case  of  A  hinges ) . 


The  recent  radial  vortices,  issued  from  blade  B,  beoome 
effective  on  blade  A  after  u  time  ZeS  ,  approximately 
equal  to  half  a  period. 


/ 
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e 


5.2412 


The  same  happens  with  the  former  radial  vortices  which 
follow  with  intervals  of  T  .  The  prooess  is  completed 
by  the  establishment  of  tlfe  final  pofeitionning  of  the 
helix  vortices  which  also  affects  blade  A  approximately 
every  ljalf  a  period  t 


During  the  visualization  tests  carried  out  for  the 
study  of  non  stationary  conditions,  very  heavy  blades 
were  intentionally  selected  (  Oj  =  0.10)  in  order  to 
make  appearant  the  aerodynamic  1  prooess  without  too 
much  blade  flapping. 

The  disc  loading  was  also  high  in  order  to  avoid  cou¬ 
pling  and  to  obtain  a  better  view  of  the  helix  vortioes, 

The  pitch  angle  disturbance  amplitude  reached,  during 
those  tests,  fifty  per  cent  of  the  Initial  pitch  angle. 

The  displacement  of  the  marginal  vortices,  after  one 
period  time,  is  not  visible  too  well  on  the  films* 

On  the  contrary,  the  displacement  of  the  blade  relati¬ 
vely  to  the  marginal  vortices  is  quite  large  if  the 
rotor  is  tilted  («*  )  .  It  is  also  noticeable  in  the 
case  of  cyclic  pitch  variation  (02);  though  at  a  smallei 
degree. 

Let  us  examine  the  case  where  the  whole  rotor  pitches 
nose-up,  by  an  angle  aL  (t),  starting  from  hovering 
conditions,  and  let  us  determine  the  response  to  thiB 
motion  in  flapping. 

The  resulting  normal  blade  section  velocity  (transla¬ 
tion)  and  the  acceleration  are  i 

Vz-  oJrcoS  V  -  r  ^  -  A  VC 

a 

Vj  s  oi  rcos'K’  -of Corbin 

The  local  lift  and  the  moment  at  the  blade  root  due  to 
the  blade  inertia  forces  and  to'  the  impulsion  of  air, 
are  : 


/. 


Aki  =o l  r  m  cos  V-  Z *5  cu  r  rr/sia  V/_  mV  d*  £ 
dr  dfc* 
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-m'rpco4*  p  12-  ('i'rcosf'-  ^  cur. tin  y_r  _  £g-AVi) 


M,  .  *  Vo5  f-i I 4,rv  V.  I.  *12.  _  I.p,  a,* 


m.2.66< 


+  f>  («d*CoJ  Y'-^cuSiaV-  ft 

(hr  *£')  -jg  +  Wf>- J  p<*wr  Tdr  + 


Av^dr 
dt  c 


[^•cos^Vp  ^3p^.aV(l#+f  /rl  avcdr 

where  m  dr  is  the  mass  of  the  blade  element  dr  and  I 
the  blade  moment  of  inertia  relative  to  the  root. 

T  +  o  7TC.*-r3 

Let  us  compare  - — —  to  G3 

0  7TCzRa  7 T  be  (T 

r  ~<?Gj  *  ~3b  “W  T* 

(  b  being  the  number  of  blades). 

l.'ith  a  solidity  (J  =  0.10  this  ratio  is  equal  to  5  %, 
in  the  case  of  a  two  blade  rotor,  that  is,  2.5  to  5  % 
of  I  for  a  conventional  rotor. 

V.'e  will,  therefore,  neglect  p  -7rC  (appearant 

mass  effect)  in  front  of  I.  “** 

’e  will  introduce  now  the  effect  of  lift  (due  to  circu¬ 
lation)  on  the  flapping  motion,  by  examining,  first, 
simple  cases,  then,  more  complex  ones. 

Usually,  the  non  stationary  conditions  and  the  free 
vortices  being  not  taken  into  account,  it  may  be 
written  : 

(III. 2. 67)  dL  s  p  uj  r  P-  _Lp  (tur)^c  dCi  ,  o<__  Co>s \y _ dfi  » 

'  4  r  dt  v  co  cor  dt  ' 


(III. 2. 68) 
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The  flapping  equation  is  then,  (f  of  (fc)  nose  up)  t 


4- to  +  2co  +■  co2,  Jfi.  b,  =  -  Ji  -  o(u>  -^S- 

2cv  *4  +  co*  <3,  -  b  -  CO  iLlrJlbj  =  -2  of  O) 

*•  Ic 


In  practice,  this  system  has  a  time  constant 
and  a  natural  frequency  to  . 


Z  I. 
(oG3 


Let  us  introduce  the  effect  of  the  recent  wake,  belon¬ 
ging  to  the  blade  from  which  it  is  issued.  This  is 
the  black-box  system  which  has  been  described  previous¬ 
ly. 

In  stationary  conditions  we  .have  : 


(III. 2. 69)  r’.itiio.fdtMl'.  £  _*£_) 

2  <U  u>  cur"  At 

The  true  circulation,  at  the  output  of  the  black-box, 
is  given  by  the  equation  which  was  indicated  in  para¬ 
graph  3 . 23  i 


(III. 2. 70) 


0.76Z 


ip-cor  T  =-  JLL  (0.639)  +°£?cv rP 


(III. 2. 71) 


For  the  calculation  of  the  flapping  motion  we  will  assum 
cor  constant  and  equal  to  the  one  at  r  =  0.75  R. 

tor  is  thus  to  be  replaced  by  W r  =  0.75  CO  R. 

Let  us  multiply  the  two  members  of  equation  III. 2. 70  by 
ptor*df  and  integrate  from  ro  to  R. 

The  flapping  equations  may  then  be  written  s 


-A 


to  r*  r  dr 


h  ’{fj  c  -ax  r,dr 


£  4  CO*  f*  .  -1L  -  £  CoS  V  -  Z  *0  jia  V 

o,639to<^  £  4  cor  co  G3  [*>  +  M  4  2^? cor  M. 

r  0.639  Co  CoS  H'+oi.tv  5ia^j4  cor  co  ci  Co*’ 


•  A. 108 

It  Is  known  that , comparatively  to  the  previous  III. 2. 68 
system, the  black-box  introduces  a  decrease  of  amplitude 
and  a  phase  shift  which  does  not  exceed  15  degrees. 

Let  us  examine  the  influence, on  the  disturbed  motion, of 
the  inner  and  outer  marginal  helix  vortices. 

Before  the  disturbance  those *helixes  are  defined  by  thelx 
circulation  intensity  and  their  pitch. 

In  the  case  of  symmetric  disturbances,  when  the  final 
conditions  are  reached,  the  helices  have  a  new  circula¬ 
tion  intensity  and  a  new  pitch. 

In  the  case  of  Antisymmetric  disturbances,  which  is  the 
case  examined  here,  the  helic  is  distorbed  in  a  sinusoi¬ 
dal  way,  its  pitch  does  not  vary  and  its  circulation 
intensity  varies  in  azimuth  (because  of  the  radial  vor¬ 
tices)  . 


During  the  transient  period  a  triple  effect  will  be  ob¬ 
tained  (in  differentiating  with  respect  to  three  varia¬ 
bles).  First  effect  :  the  marginal  helix  vortices 
remain  at  their  place  and  keep  their  pitch  and  their 
intensity.  The  blade  undergoes  a  flapping  motion  and 
takes  a  position  distant  by  Az  above  its  initial  posi 
tion.  Originally  blade  A  was  located  at  a  distance  X, 


At  a  time  t  during  the  disturbance,  the  distance  X4  of 
the  blade  to  the  first  outer  helix  is  : 

AZq,  =  K  (  (i  -  OL  COS  V  ) 


Z<L  =  *0+  »  \Jf.  X  -f-  K  C  fb-c*  CC&V  ) 

** 

and  its  distance  z i  to  the  first  inner  helix  is  * 

AZI  _  (  (i.oicosV 

2^.  z„+AZl  .  Y  +  £  (  p>  -  oi  <®S  V ) 

Let  A  V<1(  £.,(")  be  the  velocity  Induced  by  the 
helixes  in  hovering  conditions  and  £  Vi(Z0  +•  L%  ,T  )  the 
induced  velocity  when  the  blade  is  lifted  by  '  , 

We  have  : 

I  AVl  .  Az  4*;  (Z.,r) 


(z0s)  <o 

Al _ 


A 


If  Ax  >0 
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(III. 2. 72) 


Hence,  the  new  induced  velocity  variations. 

=  -2- <  ^ (i,,r) 

{  Z  AV£e.  T?  (  p-*C03  V )  M  (*  f) 

<5  V{.  -  -  (  J?>  ^  o<  cos  V)  k  -R  r  51  ViT  +•  ZA  V; 

3  y* 

The  expressions  7^  ( Z0)r )  are  derivatives  of  the 
elliptic  integral  K  and  of  the  Heuman's  A# 

function  (see  hovering  conditions). 

Second  effect.  The  intensity  of  the  marginal  vortices 
is  constant  and  their  distance  Z«*  £  to  the  blade 

remains  the  same.  The  pitch  of  the  vortices, alone, varies 
in  azimuth. 

This  term  takes  the  following  form  : 


c* 


<5vt  „  A  Yt  (HO 


W 

Third  effect.  The  pitch  and  the  distance  Z0  of  the 
vortices  are  constant.  The  circulation  varies  in  azimuth, 

This  term  may  be  written  1 


<Jvc  =  AL.  Vi  (zotr) 


The  disturbances  being  small,  we  will  neglect  the  two 
last  effects. 

Let  us  introduce  the  first  effect  (III. 2. 72)  into  the 
III. 2. 71  system. 

The  circulation  equation  is  modified  as  follows  t 


A 
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(111.3.73)  < 


0.639  G5CUG3/S+  (0.639co  02  R  k  +  0,  J62  corujGj )  p> 

+  °'l62.  cur  uj  G2  R  k  p>  +  M  + 

U  t 


0  „  .  0.^62  fju r 


M 


u  w 

=  0.639  <■<•>  Gj  «<.  cos'/  +  O.639  cu  Crj  co  sin 
+  (0^62  uyrwQj  +  0.639  ^Q2  Rk)  <i  cos  Hf 
-  0.639a>G2Rko4co  sin  ^  uJr  co  G2  R  ko(  cosK 


Introduction  of  the  term  III. 2. 72  In  the  conventional 
flapping  equations  (III. 2. 68)  may  have  a  double  effect. 


Firstly,  the  b* flapping  term  may  be  modified,  secondly, 
there  may  be  a  slight  m  spring  effect  on  the  rotor. 
This  latter  effect  improves  stability  in  hovering. 


The  effect  of  the  radial  vortices  located  beneath  the 
rotor  disc  has  been  studied  by  R.  Loewy  and  also  by 
R.  Timman  and  A.  Van  de  Vooren  (Ref.  13  and  14). 


As  far  as  there  is  no  resonnanoe  this  effect  may  be 
compared  to  a  gain  modification  of  the  lift. 


We  have  introduced  in  the  calculations  of  the  above 


paragraph  the  expressions  instead  of  the  more  accurate 


development  terms  (see  hovering  conditions).  Accoun| 
may  be  taken  of  this  by  introducing  B  R  instead  of  R 
(with  B  <  1)  in  the  higher  limit  of  the  integrals. 


3.242  -  Periodic  working  conditions 


The  circulation  P(r,t)  may  be  expressed  in  a  general 
way  as  follows  : 


m 


Ml  m 

f( r,tj  5  r  ( r,t )  *  X!  [>m  (r,t  J  cos  mH'*  r* (r,t )  41a  m  ¥ 


Only  the  first  term,  which  is  a  function  of  r  alone, 
exists  in  hovering. 


Disturbed  flight  conditions , the  knowledge  of  which  is 
necessary  for  stability  and  control  studies,  may  be  wri* 
ten  in  a  linearized  form  as  follows  : 


i  T Q  cr/k;+  <Jl}  (r,t;  cosY*  $  ^  (r,t)  sin 


The  expressions  end  TM  are,  in  periodic  working 

conditions,  independant  from  time  and  functions  ofnalonel 
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This  is  the  ideal  case  of  a  rotor  in  stabilized  level 
flight  (when  the  machine  is  stable  or  indifferent  ). 

The  phenomenon  is  periodic  :  the  marginal  helix  vortices 
have  reached  their  equilibrium  location  and  perform  only 
periodic  displacements,  when  varies.  Their  distance  z 
to  the  blade  depends  upon  the  azimuthal  position  of  the 
latter.  This  distance  z  takes  into  account  the  flapping 
effect  r  fl  .  There  is  no  reason,  therefore,  to  const- 
der  this  effect  separately,  as  this  has  been  done  in  the 
previous  paragraph  for  transient  conditions. 

The  radial  vortices,  of  the  type  studied  by  R.  Loewy,  for 
a  steady  surface  which  extends  towards  infinity.  A  theo 
retical  problem  arises  concerning  those  vortices  :  what 
happens  to  them  when  the  sheets  of  helix  vortices  coil 
up  to  form  the  marginal  voriices  ? 

From  a  practical  stand-ooint  it  may  be  noted  that  the 
intensity  of  the  radial  vortices  is  only  high  at  the 
immediate  vicinity  of  the  blade.  This  remark  permits 
approximations  to  -be  introduced. 

5.25  -  Conclusive  remarks 

The  application  of  the  non-stationary  theory  leads,  for 
fixed  wings  and  incompressible  flow,  to  satisfactory  results 
This  makes  one  believe  that  the  results . will  .also  be  good 
in  the  case  of  fairly  high  loaded  helicopter  rotors. 

In  the  author's  point  of  view,  the  primary  effect  of  the 
introduction  of  non-stationary  conditions  concerns  the  flap¬ 
ping  motion  of  the  blades  (in  transient  conditions  as  well  as 
in  the  steady  state) and  rotor  control  performance. 

The  mathematical  apparatus  which  was  introduced  in  the 
calculations  is  of  the  same  order  as  those  normally  used  to 
handle  such  problems. 


3  .3  -  Forward  flight 

(It  is  recommended  to  read  first  paragraph  3.1) 

3.31  -  Setting  in  equation  form 

Let  us  consider  a  two-blade  rotor  in  forward  flight.  Let 
1"a  be  the  non-stationary  circulation  at  a  stand-point 

( V  )  blade  A  and  1"^*  ( P,  'K )  the  stationary  circulation  a 
the  same  stand-point. 

V.'e  obtain  : 
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(III. 5 .7*0 


^V,^=^C-^(’a>r  +  V$itiV)(6c+  6fc  +  8,  cos  V'V  02  sin  </) 

"4c^rr'/siaoC* +  rTT  +  '/P><*3<0 


_  _1  c  _L  /  <?Ja  ^  _  _j_c  <*Cl  y  V. 

2  <u  4W J  Jr4  r„q  £  dc  u 

c 


The  integral  term  represents,  in  an  approximate 
flay,  the  half  helix  issued  from  blade  A. 

The  relation  between  the  stationary  circulation  and  the 
non-stationary  one,  is  given  by  : 


df7 .  2x  0.38-f  (a>e  +  \/s‘irtV) 
(III. 3 .75)  $fc+  c 


Vcv  cos  V 
CoT  +  Vsia^ 


An* 

=  0.«9  -ff  * 


’  Js^SL  ( o.r+  Vsia  V)  +0.S3S  ytuCMy1l 
l  ujr+l/sinV 


(111.3.76V 


For  -flapping*,  we  will  substitute  to  this  equation  the 
following  one  j 

+  k  T  =  0.639  kPX 

at  at 

k  =  2  *  cu  o.  75  “R 

The  flapping  equation  may  then  be  written  : 

( central,  hinge) 


(in. 3 .77)  l.p  +  cu*  1,(3  =  p  / r  (r,V  )  cvrzdr 


3.32  -  Defining  vortex  location 

As  it  has  been  said  in  the  study  of  hovering  conditions 
(paragr.3.1)  two  methods  are  proposed  : 

First  njethod. 

The  first  half-turns  of  the  helixes,  which  are  evenly 


distributed  along  the  blade  radius,  are  replaced 
by  half  straight  lines. 
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_±_  [*  ar  dfi, 

4y  Jr.  3<!  r.c 

Then  the  helix  vortices  are  grouped  together  in  four 
marginal  helixes  : 

-  two  outer  marginal  helixes  (blades  A  and  B) 

-  two  inner  marginal  heliaes  (  "  "  "  ") 

The  circulations  of  the  outer  and  inner  helixes  are 
opposed  and  equal  to  +  PM  ( <K)  which  is  the  maximum  cir¬ 
culation  to  be  determined  by  successive  approximations. 

Each  marginal  helix  is  replaced  by  a  circle  (C*)  and  a 
solenoid  (S*)  the  characteristics  of  which  are  summarized 
in  tables  3.34  and  3*35* 

TABLE  3.34 

Vortex  Circles 


A 

B 

Inner 

Outer 

Inner 

Outer 

Radius 

'•.♦ylf-'y 

R 

R 

Z 

^3(1X00+  Vlfc  (V)  T 

V',s\n.«  +  VCzi(Y)  J 

Circulation 

*rH  (r> 

+  PH  (V) 

-  rH 

Circle  center 
location 

0  B*  oosw 
• 

'0S+)T 

/ 
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Radius 


Limits  of  inte¬ 
gration 

Z  t<3  OO 


Circulation 


B 

Inner 

Outer 

H 

2 

R 

Za  -(Visiaot+VC^)  T 

+  r„ 

-r„ 

(^sirui+V^) T  (VQs\nc(^it)X 


(V0  Vtz)T 


Direction  of 
axis 


V0  sin  oLa  +  ViZi 


\A  CCS 
c 


\J0  sinU0  +  VCZt 


V0  LOS  o i 


Displacement  of 
solenoid  axis  0  A 
(see  Fig. 114) 


The  distances  it  of  the  tables  are  measured  normally  to 
the  rotor  disc,  from  blade  A  downwards. 

The  value  of  is  given  by  : 


m  *fW**A. 

The  values  of  (  Sr )  may  be  estimated  with  the  help 
of  curves  of  Figures  107  to  11J. 

The  first  value  of  r„  (V) 


(111.5.78) 


rM  +  rv  imf" 
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may  be  taken,  for  the  first  approximation,  from 
formula  III. 5. 7^  In  which  the  induced  velocity  terms 
are  replaced  by  the  meanV£trn, 

will  be  computed  for  r  -  0.9  and  for  0,  90*. 
180*  and  270*. 

It  is  important  to  take  into  account  the  circle  (C*) 
center  location.  The  solenoid  axis  displacement  (OA)  ma 
be  .neglected. 

The  effect, op  blade  A,  of  the  radial  vortices  Issued 
from  blade  B  will  be  introduced  in  the  form  of  a  correc¬ 
tive  terra,  when  the  correct  value  for  P,*  will  have  been 
obtained,  n 

The  expression  of  this  corrective  terra  is  given  in  the 
f ollowing  paragraph. 

Second  method. 

In  the  same  way  as  previously  the  first  half-turns  of 
thf  helixes  distributed  along  the  blade  radius  will  be 
replaced  by  half  straight  lines. 

However,  in  order  to  avoid  the  successive  approximation 
for  Pm  ,  the  remaining  helixes  will  no  more  be  grouped, 
here, in  marginal  helixes  but  will  be  considered  as 
sheets  of  oblique  solenoids  stretching  from  z  =  0  to 
infinity. 

The  intensity  of  every  elemental  solenoid  is  unknown  and 
may  be  written  :  ,• 


_ 1 _ 

«  r<  (  V  sirx  oi  +•  Vi.*, )  T 


The  corrective  term  for  the  radial  vortices  issued  from 
blade  B  will  be  maintained  here. 

3 .33  -  Calculation  work  up 


Let  us  write  : 


The  circulation  is  defined  by  42  coefficients. 

This  figure  may  seem  very  high,  but  it  may  be  noted 
that  this  comes  to  calculate  load  distributions  for  seven 
azimuthal  blade  positions, which  is  quite  usual  in  rotor 
calculations. 

We  may  also  note  that  flapping  calculations,  with  the 
help  of  equation  III. 3 .57,  are  simple  and  that  use  of 
equations  of  table  3.20  (still  valid)  does  not  raise  any 
difficulty,  these  equations  being  solvable  With  variable 
parameters . 

However,  the  calculations  of  the  induced  velocities  for 
variable  y  values  and  stand-points  (r)  such  as  * 


(111.3.80)  =  30°,  3 6°,  45\  60B,  72 9 ,  108°,  120°,  135",  144*  and 

150°. 

are  not  practically  feasible  without  availability  of 
Tables  of  functions. 

3.34  -  Requirement  of  Tables  of  Functions 

The  objective  is  the  calculation  of  the  Induced  veloci¬ 
ties,  for  different  azimuthal  positions, for  ¥  values 
defined  by  111.3.80  and  for  the  two  parameters  u  e  JL 
and 

toR 

When  use  is  made  of  the  first  method,  one  must  compute  : 

-  the  velocities  induced  by  an  inner  and  an  outer 
circles, 

-  the  velocities  induced  by  an  inner  and  an  outer 
solenoid. 


-  the  velocities  induced  by  the  radial  vortex  corrective 
terms. 


a. lir 


When  use  is  made  of/second  method,  one  must  compute* 

-  the  velocities  induced  by  the  solenoid  sheet, 

-  the  velocities  induced  by  the  radial  vortex  corrective 
terms. 

The  typloal  formulae  are  given  by  equations  111.3.81, 
in  Table  3.36. 

In  order  to  render  the  method  practically  accessible 
to  engineers,  calculations  should  be  carried  out  with 
the  help  of  numerical  Tables. 

Digital  computers  should  be  use<i  for  the  drawing  up 
of  such  Tables. 

3.35  -  Flexible  blades 

It  is  possible,  in  the  proposed  method,  to  take  into 
account  blade  flexibility. 

To  this  end  ,  the  value  of  the  flapping  angle  ft  must 
be  replaced,  in  equation  III.3.7A,  by  the  natural  modes  0 
a  homogeneous  blade. 

The  flapping  equation  111.3.76  is  then  replaced  by  a  linej 
system  which  will  define  the  new  7j  ^  coefficients  substi¬ 
tuted  to  the  and  b4  flapping  coefficients  1 


ccufc 


( 5ia  /Vjf 


sh/v§-; 


with  : 


J.  / 
r*  /  EI 


V  =  natural  frequency  of  the  blade 
Cl  =  flexural  rigidity. 


C/rc/tzs,  *  c/fj-Gjnc&s*  f*  with  & 
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CONCLUSION  •  • 

4.1  -  The  results  of  this  study  show  that  it  is  possible 

to  substitute  to  the  determination,  by  analysis,  of  a 
helicopter  rotor  vortex r structure,  an  experimental  method 
leading  to  a  satisfactorily  simple  and  general  vortex  ar¬ 
chitecture. 

This  experimental  method  consists  in  visualizing  the 
vortices  by  smoke  emission  on  model  rotors  operated  in 
a  wind  tunnel. 

The  accuracy  of*  this  method  is  not  very  high  but  it  pro¬ 
ved  to  be  quite  satisfactory. 

In  order  to  avoid  hasty  generalizations  of  particular  * 
cases,  a  large  number  of  operating  conditions  and  of 
blade  characteristics  was  experimented. 

However,  only  "normal”  flight  conditions,  such  as  t 

-  flight  in  established  state, 

-  disturbances  due  to  rotor  control  motions, 

-  fairly  high  disc  loadings, 
have  been  considered. 

Peculiar  cases  such  as  rotor  stall,  blade  flutter  and 
aerodynamic  resonances  have  been  intentionally  omitted, 

4.2  -  Confrontation  of  the  vortex  theory  with  test  results 
was  made  in  a  way  which  is  believed  to  be  logical. 

The  introduction  of  inner  marginal  vortices  (already 
brought  to  evidence  by  Robin  Gray)  was,  in  particular, 
considered  as  necessary  as  well  from  the  theoretical  point 
of  view,  as  for  the  obtainement  of  more  accurate  load 
distributions  over  the  blades.  , 

The  theoretical  correlation  of  the  mathematical  vortioez 
distributed  along  the  blade  to  the  marginal  vortices 
observed  in  the  win'd  tunnel  was  not  established  (same  as 
for  fixed  wings) . 

It  has  not  been  possible  to  demonstrate  with  aceuraoy 
that  the  marginal  vortices  are  to  be  considered  as  stream¬ 
lines  (in  a  co-ordinate  system  linked  to  the  blade)  and 
that  the  vortices  which  appear,  at  any  modification  of 
lift,  are  orthogonal  to  the  latter. 

Our  conclusions  refer  to  the  fact  that  the  marginal 
vortices  do  not  give  rise  to  any  exchange  of  energy 
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(  V  *  rot  V  ■  0  )  while  the  radial  vortioes  do  give  rise 
to  such  an  exchange  (1T\  rot  V  /  0) . 

In  the  Rankine  -  Froude  theory,  the  infinite  down-stream 
flow  velocity  is  twice  tHe  one  at  the  rotor  disc. 

In  the  solenoid  formulae  which  were  used, flow  velocity  has 
not  this  double  value. 

It  seems  that  a  theoretical  explanation  may  be  found  to 
this  in  the  study  performed  by  R.  Hirsoh  (see  Ref. 11). 

• 

4.3  -  The  proposed  vortex  architecture  is  intentionally  simplifie 

We  have  sought  to  avoid, wherever  possible,  giving  empiri¬ 
cal,.  and  so  called  "practical",  formulae  for  situating 
the  vortices. 

We  believe  that  such  formulae  lead  to  a  delusive  precision. 

4.4  -  Even  when  the  vortex  structure  is  available,  the  mathe¬ 

matical  problems  remain  difficult  to  solve. 

Every  time  it  was  possible  complex  functions  were  replaced 
by  approximate  expressions  comprising  polynomials  and 
trigonometric  functions. 

Such  a  "slaughter"  of  "attractive  formulae"  is,  however,  no 
always  possible. 

Different  methods  were  sometimes  suggested  for  solving 
one  same  problem.  They  are  generally  subject  to  the 
difficulty  to  handle  the  mathematical  apparatus. 

4.5  -  The  object  of  this  work  was  to  make  available  to  engineers 

a  method,  easy  to  work-up,  for  the  evaluation  of  air  loads 
acting  on  helicopter  rotor  blades. 

We  believe  that  this  result  has  been  achieved  though 
a  complementary  work  consisting  in  numerical  calculations 
remains  to  be  performed. 

The  distribution  of  the  air  loads  acting  on  a  rotor  blade, 
in  hovering  conditions,  such  as  obtained  with  the  proposed 
method  was  compared  to  the  available  J.  Rabbott's  experi¬ 
mental  results. 

It  may  be  seen  that  the  agreement  of  the  two  curves  is 
excellent,  particularly  at  the  blade  tip. 

Introduction  of  non  stationary  effects,  for  the  stability 
and  control  calculations,  leads  to  linearized  equations 
similar  to  those  normally  used  for  such  problems. 

Application  of  this  method  to  forward  flight  conditions 
necessitates  the  use  of  numerical  Functions  Tables  /• 
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which  remain  to  be  established. 

The  calculations  of  such  Tables  must  be  carried  out 
with  the  help  of  digital  computers! 

Availability  of  these  Tables  would  render  the  proposed 
method  easy  to  work-up. 

It  would  also  be  of  interest  to  havef  for  further  com¬ 
parison,  accurate  experimental  data  of  air  load  distri¬ 
butions  along  the  blades,  in  forward  flight,  with  aooount 
for  non-stationary  effects. 
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